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Abstract 

We study a generic model of self-assembling chains which can 
branch and form networks with branching points (junctions) of ar- 
bitrary functionality. The physical realizations include physical gels, 
wormlike micells, dipolar fluids and microemulsions. The model maps 
the partition function of a solution of branched, self- assembling, mutu- 
ally avoiding clusters onto that of a Heisenberg magnet in the math- 
ematical limit of zero spin components. Regarding thermodynam- 
ics properties as well as scattering strucuture factor, the mapping 
accounts rigorously for all possible cluster configuration, except for 
closed rings. The model is solved in the mean field approximation. It 
is found that despite the absence of any specific interaction between 
the chains, the presence of the junctions induces an effective attraction 
between the monomers, which in the case of three-fold junctions leads 
to a first order reentrant phase separation between a dilute phase con- 
sisting mainly of single chains, and a dense network, or two network 
phases. The model is then modified in order to predict the structural 
properties on the mean field level. Independent of the phase separa- 
tion, we predict the percolation (connectivity) transition at which an 
infinite network is formed. The percolation transition partially over- 
laps with the first-order transition. The percolation transition is a 
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continuous, non thermodynamic transition that describes a change in 
the topology of the system but not a thermodynamic phase transi- 
tion. Our treatment which predicts both the thermodynamic phase 
equiUbria as weh as the spatial correlations in the system allows us to 
treat both the phase separation and the percolation threshold within 
the same framework. The density-density correlation correlation has 
a usual Ornstein-Zernicke form at low monomer densities. At higher 
densities, a peak emerges in the structure factor, signifying an onset 
of medium-range order in the system. Implications of the results for 
different physical systems are discussed. 



1 Introduction 

Networks and branched structures are ubiquitous in both natural and syn- 
thetic materials and form under a variety of equilibrium and non-equilibrium 
conditions. In this paper we present a theory that describes in a unified way 
the structure and thermodynamic properties of equilibrium networks and 
their relation to several soft condensed matter systems, such as gels, worm- 
like micelles, microemulsions, and dipolar liquids and colloids. In all these 
systems the networks consist of cross-linked elongated objects (e.g., polymer 
chains in the case of a gel). The study of network phases is of both theoret- 
ical and practical interest. From the practical point of view, gels and sol-gel 
systems are at the core of many industrial, biological and bio-medical appli- 
cations. Examples range from novel plastics and food processing to gel chro- 
matography and tissue implants ||T0|,p6i. Despite intensive theoretical and 



experimental investigations carried out over the past three decades, many fea- 
tures of network-forming systems still remain unclear. For example, whether 
the gelation transition in the physical gels is a first-order or a continuous 
transition, is currently under debate [^,^,^,[3^. Experiments have not yet 



provided a clear-cut answer, because the gelation transition can be obscured 



by the van der Waals interaction . Another example is the transition of a 
solution of worm-like micelles to a self-assembled network that has been sug- 
gested to occur on the basis of rheological measurements; this has also been 



discussed theoretically ||5^,|53|. Network formation may also be responsible 
for the 'closed loop' phase diagrams of dipolar and magnetic liquids and 
colloids and of microemulsions E^. Cryo-electron microscopy shows clear 



evidence of coexisting network phases in dilute microemulsions [p^ . In many 
of the systems mentioned above, the chains themselves are self-assembled in 
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an equilibrium manner from a large number of monomers. 

We focus here on the thermodynamic behavior and structure of systems 
with thermoreversible cross-linking; this means that the cross-links can break 
and reform under the influence of thermal fluctuations. It turns out how- 
ever, that many of the large scale structural properties of network-forming 
phases are independent of the precise nature of the cross-links. In many of 
the complex systems mentioned above, an understanding of the local and 
large scale correlations are also of importance. Certain structural transitions 
cannot be detected by examining the thermodynamic properties alone, but 
are expressed by the nature of the correlations between different components 
of the system. One example is the percolation transition, at which a network 
spanning the whole system is formed, as schematically shown in Fig. ^ It 
is a continuous transition, unrelated to the thermodynamic properties of the 
system. Generally speaking, one would like to know both the thermody- 
namic behavior of these systems as well as such experimentally measurable 
quantities as density-density correlation functions and response functions to 
applied fields such as an electromagnetic field or hydrodynamic flow. 

To answer these questions, at least within a mean-field approach, we study 
a generic system consisting of self-assembled chains which can branch and 
form networks. Each chain consists of a large number of 'monomers'. The 
chains are self- and mutually- avoiding but are permitted to branch (cross- 
link). We shall call both the branching points and the cross-links 'junctions'. 
The exact physical interpretation of the 'monomers' and 'junctions' will dif- 
fer form system to system. In each particular case, the energy and other 
properties of the junctions and the ends can be calculated from microscopic 
considerations, e.g., molecular packing for surfactant systems, or dipolar in- 
teraction energies applicable to ferrofluids. 

We present here a lattice model which establishes the equivalence between 
a solution of branched self-avoiding self-assembling chains and a Heisenberg 
magnet in the limit of 'zero' spin components; we show that their partition 
functions can be mapped onto another. This model, known as 'n = 0' model, 
was proposed by De Gennes to study polymer solutions and we have modi- 
fied it in order to include the possibility of branching. Regarding the ther- 
modynamic properties, including the density fluctuations, the model is exact 
and accounts for all possible configurations of equilibrium branched clusters 
(excluding separate closed rings, whose influence is or importance only at 
very low densities). This formulation of the problem enables us to explore 
both thermodynamic and structural properties of self-assembling branching 
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chains in a unified manner. The model has been treated in the mean field 
approximation and our major conclusions are summarized in Fig. |I|, which 
shows the phase diagram of the system as a function of the density and 
temperature. We find that despite the absence of any specific interactions 
between the monomers, the presence of the junctions induces an effective 
attraction between the monomers. The model can be amended in order to 
study the topological properties of the system and correlate them with the 
thermodynamic behavior. As the monomer density is increased or the tem- 
perature decreased, the system passes through a percolation (connectivity) 
threshold, where a network spanning the entire volume is formed. This tran- 
sition is purely topological, and has no thermodynamic signature. However, 
the junction-induced attraction does modify the concentration of monomers 
at which the percolation threshold occurs; as the temperature is decreased, 
the threshold is decreased compared with its value in the limit of infinite 
temperature, where the interactions are irrelevant. 

For three-fold junctions, the junction-induced attraction is strong enough 
to drive a first-order phase separation, where the system separates into a 
low-density and a high density phases. This transition is of purely entropic 
origin, because there are no any specific interactions between the monomers. 
The physical reason for the transition is the higher entropy of the junctions- 
dominated high density phase: although the translational entropy of the 
chains is lower in the high density phase, it is overcompensated by the in- 
crease of the entropy of the self-assembled junctions, abundant in the high- 
density phase. In this respect, this entropy-induced transition is similar 
to the crystallization of rigid spheres and the isotropic-nematic transition. 
The transition line terminates at a critical point. For four-fold and higher 
functionality junctions, the junction induced transition is too weak to drive 
a phase separation. In this case, the junction- induced interaction merely 
renormalizes the excluded volume interaction between the chains, and drives 
the system closer to the G-point. The structural percolation transition, of 
course, is also present above the critical temperature of the first-order phase 
separation. At very low temperatures and monomer densities, the chains are 
depolymerized and the system consists mainly of separate monomers. 

The nature of the phases in equilibrium depends on the temperature, but 
also on the rigidity of the chains and details of the junction configurations. 
In general, there are three possibilities: (i) a phase of dilute chains that 
coexists with a connected network, (ii) two coexisting networks, or (Hi) the 
coexistence of dilute and dense phases of disjointed, branched aggregates, 
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although not predicted by our model, is also possible. 

Treatment of the spatial density variations within the mean field theory 
shows that the density-density correlation function, relevant for scattering 
experiments, has a simple Ornstein-Zernicke form for relatively low densi- 
ties of the self-assembling monomers. For higher densities, a peak emerges 
in the structure factor as a function of the scattering wavevector; this indi- 
cates the emergence of medium range (longer range than the lattice size but 
not long range order) correlations in the system. In the intermediate regime 
the structure factor is a monotonically increasing function of the wavevec- 
tor. This is expected in any dense system: at high densities, the system 
has a low compressibility at long wavelengths and this suppresses the small 
wavevector scattering. The absence of a peak in the structure factor does 
not imply the absence of the network: there is a region in the phase diagram 
beyond the percolation line ( where an infinite network exists) where there 
is indeed no peak in the density- density correlation function. The precise 
location and strength of the peak depends on the monomer density and the 
number of junctions present. As discussed in Section. |^, these theoretical 
predictions may be related to the scattering peak observed experimentally in 
bicontinuous microemulsions. 

The predictions of the thermodynamic and structural properties can be 
put in the context of previous work on the physics of 'living' (self-assembling) 
polymers and networks; this area has been the focus of extensive experimen- 
tal and theoretical attention over the past two decades. Most notably, phe- 
nomenological Flory type (chains with linkers) theories have been employed 
to model physical gels Another important contribution is the work 



of Drye and Gates that was motivated by studies of worm like micelles BT 



However, the predictions of these various studies regarding the properties of 
the gelation transition vary depending on the details of the models that are 
studied, the assumptions regarding the solvent quality and other parameters. 
This is partially due to the fact that these studies [|29|,^ did not emphasize 
the extreme sensitivity of the gelation transition to the functionality of the 
junctions and the importance of the competition of the branching points and 
the free ends. These previous, more heuristic, mean-field theories could not 
be extended in a simple manner to predict the correlations. 

We have therefore focused on the more rigorous 'ra = 0' model, which 
allows one to evaluate exactly the number of all possible configurations of 
branched self- and mutually avoiding chains on a lattice. Lattice approxi- 
mation models the configurations of the chains only approximately at short 
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lenghtscales. However, as long as one is interested the behavior of the system 
on lengthscales much larger than the lattice constant, (i.e., thermodynamics) 
the choice of the lattice does not influence any of the qualitative predictions of 
the model. The 'tt, = 0' model has been successfully used to explain the prop- 
erties of associating monomers and has been applied to the polymerization of 
sulfur and worm like micellar solutions . The extension of the model 
to n = 1 was used to model systems consisting of chains and rings in equilib- 
rium Our model is similar to the one employed by Isaacson an Lubensky 
p3| to which it can be related (cf. Appendix D). These authors studied the 
continuum version of the 'n = 0' model with a cubic term in the Hamiltonian 
in the context of percolation focusing on the scaling behavior of the model to 
study the properties of the gelation-percolation transition but did not calcu- 
late the thermodynamic transitions (such as the phase separation predicted 
here) or study the monomer density fluctuations and scattering structure 
factor. In our model, however, the junctions are not modelled as point-like 
objects, and the treatment of topological structure is simpler than of Ref. 
||23|| , but is applicable only in the mean field approximation. The n = 1 
model, which accounts also for closed rings, but does not properly take into 
account the self-avoidance, was used by Lequeux, Elleuch and Pfeuty to 
study general properties of micellar solutions without considering the topol- 



ogy of the network. Similarly, in the two-field model of Ref. PT| , correlations 



and structure were not studied and no analytic expression for the free energy 
was obtained. 



This paper is organized as follows. In Sec. |2.1| we review the 'n = 0' model 
without junctions. In Sec. ^]2| we extend the model to include the possibility 
of junctions and establish the correspondence between this model and the 
physics of the problem of self-assembled networks. In Sec. ^ we calculate the 
free energy and phase diagram in the mean-field approximation. We next 
extend the mean field theory to include spatial variations of the monomer 
density and in Sec.^,predict the spatial variations and correlations. In the 
Sec. ^ we show how the model can be modified on the mean-field level 
to enable us to predict the topological properties and correlate them with 
the thermodynamic behavior. The approximation employed, is of the same 
level as the classical Flory-Stockmayer theories of gelation. In particular, in 
Sec. we discuss the structure and evolution of a single branched cluster 
as a function of monomer density and the temperature, and the emergence 
of a connected network. Sec. ^ describes the extension of our model and 
its predictions to the case of rigid and semifiexible chains. In Sec. ^ we 
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Figure 1: Phase diagram of solution of self- assembling branched chains as a 
function of temperature T in units of the end energy eg, and monomer volume 
frcation cf). The junction energy has been chosen to be = ee/4. The thick 
line shows the phase separation region. The dashed line is the percolation 
line. To the left of this line the system consisits a solution of disconnected 
(but possibly entangled) branched chains; to the right there is an infinite con- 
nected network. Moving along the line 6, the actual gelation transition will 
be seen as a continuous non-thermodynamic transition, while moving along 
a will result in the first-order thermodynamic transition. Thus, coexistence 
of two connected networks is possible at high enough temperatures. 



discuss the results and their application to particular physical systems such 
as physical gels, microemulsions, dipolar fluids and wormlike micelles. 



2 Zero-component' Heisenberg model 
2.1 Self-avoiding chains. 

For systems with no junctions, the 'n = 0' model has been extensively studied 
in the context of polymer solutions and micellar systems P,|l^,[l^,|T^. For 
the sake of completeness we review here the derivation of the model; in 



Sec. p.2| we extend the model to include junctions. Readers familiar with 
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'n = 0' model can proceed directly to Sec. |2.2|. We consider a system that 



comprises a solvent and a collection of self-assembling chains, each consisting 
of a large number of identical monomer units. The physical nature of the 
'monomers' can be different. For example, in physical gels these would be 
the individual molecular units that comprise the polymers, in microemulsions 
these would be surfactant covered oil (or water) domains and in dipolar 
colloids, the colloidal particles (cf. the discussion in Sec. |^). The chains 
are self- and mutually avoiding, that is, they cannot intersect themselves 
and each other. We focus on self-assembling systems where the monomers 
are relatively weakly associated within a chain; thus, individual monomers 
can be freely exchanged between the chains. In this case, the equilibrium 
distribution of chain lengths is not fixed, but is determined by the values of 
relevant physical parameters, as we discuss later. We will consider the system 
in the grand-canonical ensemble, in equilibrium with the bath of monomers 
of the chemical potential /x. We assume that the energy of a monomer at 
the end of the chain, eg, relative to that of a monomer in the middle of the 
chain is positive, > 0. This is indeed the case for mi cellar systems and 
microemulsions, because (in the regime where long chains are formed) the 
bending energy of the surfactant layer is lower in the cylindrical part than 
in the semi-spherical end-caps [^. The same is true for dipolar particles, 
because the electrostatic energy of a dipole in the middle of the chain with 
two neighbors, is lower than that of a dipole at the chain end, where it has 
only one neighbor dipole. If eg < 0, the self-assembling chains will be very 
short and the solution will mostly consist of individual monomers; this is not 
the limit that we are interested in. The end energy, taken with the opposite 
sign, -eg, can also be interpreted as /i^, the chemical potential conjugate to the 
number of chains in the solution, Np. This interpretation is usually adopted 
in discussion of regular polymer solutions 0,0. For the systems of interest, 
the energy eg is determined from microscopic considerations, be it molecular 
packing or electrostatics. The number of chains and average molecular weight 
are not fixed but determined by physically controllable parameters, namely, 
the temperature, the total monomer density (p and eg. For simplicity we 
imagine that the monomers and the solvent molecules occupy the sites of a 
(i-dimensional lattice. With the above notation, the grand canonical partition 
function of a solution of equilibrium chains is given by 

Z= J2 exp[f,N/T]exp[-eMT]Af{N,N,) (1) 

{N,Ne} 
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where N is the number of monomers in a given reahzation, is the total 
number of free ends, and Af{N,Ne) is the number of ways to arrange Ne/2 
chains of total length on a lattice. Thermodynamic quantities can be 
obtained from Eq.(|l]) by differentiation with respect to the relevant param- 
eter. For example, the mean number of monomers and chains, N and Np 
respectively, are: 



N 



d\nZ 
d\nZ 



The direct combinatorial calculation of Af{N, Nf.) is a formidable task. 
A way of circumventing the difficulties involved in its calculation, was pro- 



posed by De Gennes |T^. It was later refined by equivalent methods by 
different authors The calculation of ^/{NjNe) is done by relating 

the partition function of the Eq.(|I]) to that of a certain spin model on a lat- 
tice. Consider the ra-component Heisenberg model on a ci-dimensional lattice 
whose Hamiltonian is given by 



H 



^ ] J Si ■ Sj ^ ] hi ■ Si 

{s,} {s,} 



(2) 



where the summation is over all distinct pairs of nearest neighbors i,j. Spin 
Si is an n-component vector normalized at each site as 



n 



where a indexes the components of Si. We arbitrarily choose the field h to 
point in the direction (1, 0, 0). 

The partition function corresponding to this spin Hamiltonian is 



Zh = Tr exp[-H/T* 



(3) 



where the trace operator Tr signifies an integration over all possible direc- 
tions of Si, divided by a normalization, so that the trace Tr A of any quantity 
A is: 



Tr A 



I n. da 



(n) 
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where ^2*^"^ is the n-dimensional sohd angle. The 'temperature' T* has no 
physical meaning and is not related to the real physical temperature of the 
system of in Eq.(|I]); we put T* = 1 (another choice of T* would amount to 
a renormalization of the constants J and h). Expanding the exponential in 
Zh of Eq.(|D in a power series, we obtain 

= E n('^^^^.- + (4) 

Each term in Eq.(^) has the following form 

Tr \^J^ SiS j ...Sii S jiK"" Sp^\...Spi (5) 

m times n times 

Note that in the tracing operation in Eq.(^, all sites are decoupled. 

In any realistic system in which the spin S with n components is to 
have a physical interpretation, n > 1 is required; n = 1 corresponds to the 
usual Ising model. However, quantities such as partition function and the 
cumulant expansion remain mathematically meaningful even when n < 1. In 
particular, one can consider the mathematical limit n — > It can be 

shown that when n — > (cf. Appendix |8.1|) , 



Tr Si^aSj^p — Sa/sSij (6) 

and all other cumulants are zero. The spin S does not have a physical inter- 
pretation in the n —>■ limit and should be regarded purely as a convenient 
mathematical device. Precise mathematical meaning can be given to Eq.(P) 
(cf. Appendix |8.1| ). In particular, Eq.(^) means that in the expansion (^) the 
only non- vanishing terms are those in which every spin Si appears twice. As 
can be seen from Eq. (^) , this condition is satisfied by terms consisting of the 
chains of neighboring bonds: ...JSiSjJSjSj'.... Since every spin must appear 
twice for the term to be non-zero, the uncoupled spins must be paired off 
with another spin at the same site either by (i) closing a chain on itself or 
by (a) the appearance of the single-spin, 'field' term hiSi. However, in the 
limit of — i> 0, the terms given by (i) vanish, due to summation over all 
the components (via summation over the index a); this sum is proportional 
to the number of components, n, which is equal to zero. The terms that 
contain the field variables, hi, do not vanish because the /i-terms single out 
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one component a, parallel to the field h. It is easy to see that such terms 
consist of products of the form 

Tr hi^iSi^iJSi^iSj^iJSj^iSk^i---JSq^iSr^ihjiSr^i (7) 

that repeat themselves with different i and r. Any term of the form of Eq.(^ 
corresponds to a self-avoiding random walk on a lattice, starting at site i and 
ending at site r. In other words, a term containing the factor J"^h?^ counts 
all the possible configurations of k random walks of total bond length m. 
Thus, the partition function can be written: 

Zh= J'^'h^'J^iN^.N,) (8) 

where N'{N, Ne) is the number of ways to arrange on a lattice, an ensemble 
of self-avoiding random walks with a total number of bonds Nb and a total 
number of ends Ng. Noting that the number of chains is Np = Ne/2 and the 
number of monomers is N = Nb + Np, one can see that Eq.(|D is identical to 
Eq.(|l]) if the following identification is made 

.1 = 6" 

h = e-'^J'/^ = hoJ'/^ 

Each term in the sum in Eq.(^ corresponds to a different realization of the 
grand canonical ensemble of a solution of self-avoiding chains. We have thus 
demonstrated that the grand canonical partition function of a solution of 
polydisperse, living polymers is identical to that of the n-component Heisen- 
berg model where the number of components — > 0. As mentioned in the 
Introduction, the lattice approximation is adequate as long as the chains are 
much longer than the lattice constant. This theory describes the ensemble 
of self- assembling 'living' chains, which we are of interest to us. The molec- 
ular weight and the length distribution are not fixed but are a functions 
of external parameters, such as the monomer density. However, this model 
has also been successfully applied to non-self-assembling polymer solutions 
(where the degree of polymerization is constant and fixed by the chemical 
preparation technique). The continuum version of this model is the basis for 
the application of field-theoretical methods to polymer physics |T^,[T^,|5S 



The self-assembling theory can be used to predict the properties of non-self- 
assembling polymer solutions of long chains (small h in our formulation), 
because the scaling behavior of such systems is universal and independent of 
the detailed chain properties. 
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2.2 Three-fold junctions 

The model presented in the previous section can be modified to include the 
possibility of junctions that connect several chains; this allows a formation 
of a branched structure and is achieved by introducing an additional term 
to the Hamiltonian of Eq.(^. Again, we consider the grand canonical par- 
tition function of an equilibrium solution of branched cross-linked polymers 
with reversible cross-links. The cross-links (branching points) are in thermal 
equilibrium and can break and reform. Similar to the solution of polymer 
chains discussed above, the grand canonical partition function of the system 
is 

= Yl exp[^iV/T] exp [-eeNjT] exp[-e,iV,/T] Ar{N, iV„ Nj) (9) 

{N,N,,Nj} 

where is the energy of a cross-link relative to the energy of the monomer 
in the middle of a chain, Nj is the number of junctions (cross-links) and 
N{N, Ne, Nj) is the number of ways to put self-avoiding branched chains of 
total length on a lattice so that Ng ends and Nj junctions are formed. The 
chemical potential for the monomers, /i and the end energy eg are defined as 
in the previous section. Similar to the discussion of the end energy and its 
associated chemical potential, the junction energy taken with the opposite 
sign, -ej, can be interpreted as the chemical potential conjugate to the number 
of junctions. We suppose that the physically controllable parameters are the 
monomer density 0, the temperature T and the ends/junctions energies ej. 
In this formulation, the density of ends 0g and the density of junctions (pj 
are not fixed but determined in thermal equilibrium as functions of (p and T. 

In this section we show how the three-fold junctions can be described 
by 'n = 0' model. The generalization to junctions of arbitrary functionality 
is discussed in the next section. As in the previous section, one can relate 
the partition function of Eq.(^ to an equivalent spin model on a lattice. 
The mapping accounts rigorously for all possible configurations of branched 
self-avoiding clusters, except for closed rings. The infiuence of rings on the 
properties of living polymers has been studied in Ref. [^]. However, in 
the case of branched chains, the infiuence of closed rings on thermodynamic 
properties of the system is small, except at very low densities, and they 
are not treated in this paper. The reason for this is that for any closed 
ring there is an exponentially big number of branched clusters formed by 
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attaching side chains to it. No topological information regarding, e.g., the 
network formation, or the cluster size distribution can be extracted from 
the model formulated below. However, it can be modified to include the 
treatment of topological properties too, as shown in Sec.|^. 

We introduce the following Hamiltonian that contains a term which cou- 
ples three adjacent spins: 



H, = J2h-S, + J2JiSi- Sj) + J2 KSi,iS^,,S^,k (10) 

i ij ijk 

where i, j, k sum is over all distinct triplets on the lattice. The field h is chosen 
to point in the direction (1,0...0). The partition function corresponding to 
this Hamiltonian is 



ZH,=TT{s.}exp[-Hs{S,}/T*] (11) 

The 'temperature' T* can be taken equal to unity and any other choice would 
amount to a redefinition of the constants J, K and h. Proceeding as in the 
previous section, we expand Zh^ in series in powers of K, J and h. In the 
limit n —>■ the terms in the expansion of Zhj in powers of J, K and h 
are similar to those discussed in the previous section. However, there are 
additional terms obtained by inserting into any term of the expansion of 
Eq. (^ combinations of the form 



.JSj'^lSj^l \ KSj^lS m,lSn,l ]jSn,lSn',lJSm,lSm',l--- (12) 



junction 



Each insertion of this kind corresponds to a three-fold junction on a lattice, 
joining the sites j,m,n (cf. Fig. 0). Thus, a general term in the expansion 
of Zh-j might look like 

Tr hSi^iJSi,iSj,i[KSj^iSn,lSk,l]JSn.lSn' .ihSn\i---JSk^iSm.ihSm,i (13) 

' ^ ^ ' 

junction 

This particular term, for example, corresponds to a three-fold branched, self- 
avoiding chain, which has a junction between the points j, n, k, and free ends 
at points i, n' and m. In other words, the partition function is 
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Zu, = Tr exp[-Hs] = Yl J'^'h^'K''^ XiNb,Ne,Nj) (14) 

{NbN,,Nj} 

where M'lNh, N^, Nj) is the number of ways to arrange on a lattice an ensem- 
ble of self-avoiding branched random walks with a total number of bonds Nb 
(excluding three 'ghost' bonds at each junctions, cf. Fig. |]), a total number 
of ends Ne, and a total number of junctions Nj. The closed rings fall out 
form the expansion |l^ due to summation over spin components. Each ring 
produces a contribution proportional to n which tends to zero as n —>■ 0, 
analogously to the case of linear chains. From simple geometric considera- 
tions, the number of monomers, N, is always N = Nf, + ^N^, + ^Nj (we do 
not count the junctions as monomers; see Appendix \i.'2\ ), regardless of the 
presence or absence of loops. Note that this formulation does not treat the 
junctions as point-like objects. Choosing different forms of three-spin cou- 
pling in Eq.(|10]) one can study the influence of the junction size on the system 
properties (cf. Appendix |8.2| ). If the following identifications are made 



J = e^l^ (15) 

K = e-'^/^ji = KqJ^ 

then the Zh^ of Eq.(|l^ becomes identical to Z^ of Eq.(P). To summarize, 
we have shown that in the n — limit, the grand-canonical ensemble of 
a solution of branched, reversibly cross-linked polymers is equivalent to a 
Heisenberg magnet with three-spin term. It is important to emphasize that 
although the closed rings ( linear, unbranched, chains closed to form a ring) 
are not included in the expansion ([T3|), all the intracluster loops are counted 
properly. 

The concentrations of the monomers, ends and junctions, 0, 0g and (pj 
respectively can be obtained by differentiating the partition function with 
respect to a relevant parameters as follows from Eqs.(ptp!4D: 



1 d\nZ , 1 dAnZ 



V d\nJ' VdslnJ 
1 d\nZ , 1 dlnZ 



(16) 



V d\nh' VdlnK 



14 



where ds denotes a derivative with respect to J with the generahzed fugaci- 
ties K, J and h taken to be independent. The great advantage of the present 
formulation is that it allows to apply the methods of statistical mechanics 
developed in the past four decades for the treatment of the spin models. 

In some systems, junctions may form spontaneously; in that case, the 
number of junctions, distribution of the chain intervals between the junc- 
tions and the branched clusters' size will be determined by the junction 
energy ej, the monomer density and the temperature T. This is what one 
expects, for example, in dipolar fluids or microemulsions. In other systems, 
junction formation may only be possible if linker molecules, that connect 
several chains, are added to the system. In this case, the grand canonical 
formulation presented here amounts to the assumption that the system is in 
equilibrium with a bath of linker molecules with chemical potential -ej. 

The grand canonical potential per unit volume, uj{ej,ee, fJ^), is given by 

uj{ej, ee, fi)/T = In Z{Ko, ho, J) (17) 

Other physically realizable situations are related to Eq.(^) by a Legendre 
transform. For instance, the free energy corresponding to the case where 
junctions form only in the presence of linker molecules is 

g{fj., ee, Nj) = u - In Kq 

Of course, in the thermodynamic limit {V — > oo) the properties of the system 
are identical in either ensemble. Here, we will focus on the case where the 
physically controllable parameters are the total monomer density, 0, and the 
defect energies, eg/T = — ln(/io) and ej/T = — In(i^'o)- The corresponding 
Helmholtz free energy per unit volume is 

f = uj + fi(j) = uj + (f)ln J 



3 Mean field theory 

3.1 Mean field approximation 

The first approximation in evaluating the free energy is to disregard the spa- 
tial variations of the densities and the long range correlations between the 
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Figure 2: Junction can be formed between three adjacent sites i,j,k. The 
Boltzmann factor for a junction is Kq = e~'^^^'^. 

monomers. In this section we assume that the system is spatially uniform 
and calculate the free energy as a function of the average monomer volume 
fraction (p, and the ends and junction densities, 0g and (pj, respectively. A 
mean-field calculation of the free energy for a system with no junctions has 



been presented in Refs. p9|,p0|. We extend this calculation by an alternative 
method to a system that contains both ends and junctions. We first calculate 
the grand canonical potential f2/T = — In Z in the mean field approximation, 
which is generally known to predict the qualitative thermodynamic behav- 
ior correctlyip. As we shall show subsequently in Sec. |^, the mean field 
approximation neglects the internal loops in the branched clusters, because 
only local properties are preserved in the mean field approximation. Ge- 
ometrically, it is equivalent to Flory construction [^] when each cluster is 
constructed by adding consecutive bonds, disregarding the positions of the 
previously placed ones, and the long range correlations are lost. In particu- 
lar, intra-cluster self-avoidance is neglected in the mean field approximation. 
However, the excluded volume between different clusters is taken into ac- 
count. The partition function Z is given by 

Z = Tr exp [ - ^ JSi -Sj-J^ ^^i.i^i,,^!,, -Y.^-S] (18) 

i,j i,j,k i 
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The sums in the argument of the exponent of Eq. (|T^ are over distinct pairs 
and triples, as explained in Sec. |2.2| . We now note the following identity 

SiSj = {Si-S){S,-S) - 52 + S{Si + Sj) ~ -S^ + S{Si + Sj) 

where S = (Si). The underlined term in the equation is quadratic in the 
deviation of the local spin from its average value S. It is this term that is 
neglected in the mean field approximation. As the 'magnetic field' h has been 
chosen to point in the direction (1, 0..0) the only non-zero component of the 
average spin is Si^i. Therefore, the transverse components of the spin S do not 
contribute to the mean field approximation, due to the fact that {Si^±) = 0. 
So, in the following we drop the vector sign and Si signifies the component 
of Si in the direction of h. Similarly, the mean field approximation for the 
three-spin term is: 



s,s,Sk ~ -253 + S{S, + S, + Sk) 

Thus, the mean field approximation to the Hamiltonian reads 

Hmf = V{-qJS^ + 2aKS^) - J^i^'^^ + ^(^KS^ + h)Si 
2 

i 

where V is the number of lattice sites (which tends to infinity in the ther- 
modynamic limit) and the sum is over all lattice sites; q is the lattice coor- 
dination number. The prefactor a = where a' is the number of possible 
configurations of j, k in the triple term i,j, k making a junction around the 
cite i. The exact value of a' depends on the type of the lattice used and on 
the physical assumptions made in counting the physically relevant junction 
configurations. Physically, it represents the entropy associated with local 
rearrangements of a single junction and we discuss its value for semi-flexible 
or rigid systems in Sec.||. In most of this paper, we shall use the value of 
a' appropriate for a simple cubic lattice a' = q{q — 2)/3 (cf. Appendix p.2| ). 
A particular feature of the n — limit, following from Eq.(|D [jl], |12|, |l9l (cf . 



also Appendix 8.1) is the fact that for an arbitrary vector k, 



limTr e^^' = 1 + -^ 

n^O 2 
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Consequently, performing the trace in Eg. ([TB|) , we get 

^{n/T) =uj = ^qJS^ + 2aKS^ - In (1 + ^{qJS + 3aKS^ + hf) (19) 
The average value of the spin, S, is given by 

i 

The total monomer density is 

13 dsUJ 1 du 3 duj duj 



as follows from Eq.([16|), and can be also understood by a simple geometrical 
argument (cf. Appendix p.2|) . Finally, combining all the equations we find: 



u = ^qJS^ + 2aKS^-\n{l + ^{qJS + 3aKS^ + h)^) 
duj qJS + 3aKS^ + h 



dh l + l{qJS + 3aKS'^ + hy 
- S { qJS + 3aKS^ + h) 



Denoting qJS + 3aKS^ + h hj x, the solution of this system of equations 
gives 

(h = -Sx 
^ 2 

S - ^ 



1 + (j)x/S 



Consequently 



x 



oj = gJ0(l-0) + 2air(20)^/2(l- 0)3/2 + ln(l 
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It follows that in the limit of systems with a relatively small number of 
junctions and ends, compared to the total number of monomers in the chains, 
(/io(20)V2«0, iro(20p«0), 

= (a; _ /, _ ^aKS'')/S ^ -^(1 - - SaJ^o v^) (21) 

1 — © \/2o 
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Because we are interested in those systems where the density, 0, (and not the 
Lagrange multiplier, J, ) is a physically controllable parameter, we perform 
a Legendre transform / = + In g J in order to obtain the Helmholtz free 
energy per unit volume /(0, /iq, -K'o)- Up to the first order in ho and Kq^wq 
find: 



//r.,i-,)(>„a-,)-i)-^^-«! 

where q is the coordination number of the lattice. The first term in Eq.(p2D 
describes the self-avoidance of the chains. The second and the third term 
are the densities of ends and junctions, respectively, as can be seen from 
Eqs. ([T6|) . Denoting the volume fraction of ends by 0g and of junctions by 0^ 
it follows: 

^o(20)^/^ , 

qi/2 = '^e (23) 

g3/2 

The expansion in powers of ho and Kq is justified because we are inter- 
ested in the limit where the density of ends is much smaller than the total 
monomer density, 0g -C 0. The density of junctions is always smaller than 
the total monomer density, 0^ <^ as follows form Eq.(^) due to the fact 
that Ko = e-^^l^ < 1 for > 0. 

We note in passing that the choice of the lattice affects only the numerical 
prefactors in Eqs.(PB|),(P^). The lattice-specific dependence of junctions and 
ends densities in Eq.(^) on q is consistent with a simple probabilistic ar- 
gument. In equihbrium the probabilities of defect (e.g., ends and junctions) 
formation and break-up must be equal. The probability of bond formation is 
proportional to ^(p^q, because a bond can be formed whenever two ends are 
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neighbors on a lattice (the factor | is due to indistinguishabihty of any two 
ends). The probabihty of bond break-up is proportional to the total number 
of bonds A^;, ~ A^. Taking into account that the formation of two ends from 
a single bond costs an energy 2ee, the probability of a bond break-up is 
g-2ee/T Equating the probabilities of bond breaking and bond formation pro- 
duces Eg . (|23|) (recalling that ho = e"^*"/^). Similarly, three ends can coalesce 
to form a junction, giving in equilibrium (pl/e'^^" ^ a(f)j/e~^^ ( Kq = e'^-"'^). 
Note that this argument is independent of the fact whether the 'ground state' 
consists of infinite chains or closed rings of linear chains. 

Together with Eg. (|23|) , the free energy / of Eg.(p2D can be cast into the 
lattice-independent form 

//T = + (1 - 0) ln(l - 0) - 0, - 0,. (24) 

The first term in Eg.(p^) is due to self-avoidance between the chains. This 
particular form of / is expected on very general grounds. If one thinks 
of the ends and junctions as defects in the system of infinite chains, each 
defect lowers the free energy by ksT, which is true for any system with 
non-conserved defects, e.g. dislocations in crystal structures. 

3.1.1 Average distance between the defects 

The average volume fractions of junctions, 0^, and ends, 0g, are not fixed 
but depend on the total volume fraction of monomers, 0. Similarly, the 
self-assembling nature of the system means that the size distribution of the 
branched aggregates is polydisperse. 

Absorbing the lattice-dependent prefactors in Eg.(p3D into the definition 
of the constants K and h, it is rewritten as: 

(f). = K'off)^^ and 0, = /i'o0^/^ (25) 

As can be seen from Eg. (|25|) for <^ hQ/KQ,the number of junctions is 
much smaller than the number of ends, 0^ <^ 0g while in the opposite limit 
0j ^ 0e- Let us consider the mean length of a chain segment in between two 
consecutive defects (i.e., ends or junctions), I. Each junction is attached to 
three chain segments while each end is attached to a single segment. The 
total length of the segments is proportional to the total volume of monomers 
in the system. Conseguently (|Aj- + ^N^)! = N; the factor | corrects for the 
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double counting of each segment. Consequently, we find: 



- 20 20 

30,. + 0, 3/^^(20)3/2 + /,^(20)V2 I ) 

For -C Hq/Kq, the ends dominate and the mean chain length in between 
end points is T ~ 0^^^ while for ^ h'^/ K'^ the junctions dominate and the 
mean chain length in between junction points is T ~ (f)~^^^. These results 
are in agreement with the known results for both non-branched micelles and 
pure networks These results suggest, that at high densities (where 

the ends are negligible), a connected network is formed, as we shall actually 
prove in Sec. p.2| . They also imply that the mean chain length between 
defects, /(0) has a maximum around ~ h'^/ K'q. One can easily convince 
oneself taking the derivative ^ that I is indeed a non-monotonic function of 

1 h' 

whose maximum is located at = i^^w other words, where 0g = 30,. 
This fact will become important for the determination of the details of the 
evolution of the system from a state of disconnected chains to that of a 
connected network (cf. Sec. |5^ ). 

In the absence of junctions (e.g. when the junction energy is very high and 
the resulting value of K is close to zero), the mean chain length in between 
end points is iV = = /iq0"'^/^. The free energy F(0, ho) of Eq.(P^) can be 
transformed to a free energy that depends only on and 0g (or the mean 
chain length N) by means of the Legendre transformation /(0, /io)+0e In /iq = 
f{(j),N). The result is identical to the Flory-Huggins expression for the free 
energy of the polymer solutions: 

/(0,iV) = (l-0) ln(l-0) + iln0 

As mentioned above, this is because the mean field approximation is geomet- 
rically equivalent to Flory lattice construction, which disregards the intra- 
chain self-avoidance, but retains the excluded volume interactions between 
different chains. Thus our model describes the case of polydisperse chains 
with both ends and junctions as well as the limiting case of finite chains (with 
no junction points) where the average length is well defined; this is the case 
that is applicable to chemically prepared polymeric chains. 



3.2 The junctions-ends transition 
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The free energy / of Eq.p^) is unstable in a certain range of parameters ( 
monomer density 0, defect energies e^, and and temperature) and shows a 
first-order phase transition, terminating at a critical point. The spinodal line 
associated with this transition is determined by the condition = and 

at the critical point, = 0. We determine the conditions for the critical 
point from our model as: 



-3/2 



Recalling that ho = e~^''^'^ and Kq = e"*^^/^, the Eq.(p7D can be solved 
analytically for small 0, expanding (1 — 0) (In (1 — 0) — 1) |(/)^, which 
gives 

T. = r^^^T^ (28) 



In q^/{4:a 
-e 



1 •? „(€j-ee)/T, 



2a 

while for arbitrary they can be solved numerically (Fig.|^). It is interesting 
to note that in the small density limit, 0g ~ 0^ at the critical point, which 
emphasizes the fact that the transition is junctions-induced. For T > the 
system is a homogeneous mixture of chains and branched aggregates, while 
for T < Tc there is a two-phase equilibria between an end-rich phase that 
coexists with a junction-rich phase. As will be shown in Sec. p.2| , the end-rich 
phase usually consists of dilute, disconnected chains, while the junction-rich 
phase is usually a connected network that spans the whole system volume. 

Although there are no direct interactions between the monomers in our 
model, junction formation induces effective attraction between the monomers, 
and it is this effective interaction that leads to the phase separation of 
Eq.(P^ . The Eq.(^) has several interesting consequences: First, the phase 
separation is possible only when ej < |ee (as long as g^/(4Q;^) > 1); in the 
opposite case, the junctions are present only as a minor 'perturbation' to 
a system of linear chains and their number is too small to generate an at- 
traction large enough to drive a macroscopic phase separation. Second, the 
critical monomer density 0^ < |; this fact is important for the understanding 
of the correlations between monomers and between junctions in the system. 
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Third, if a is large (g''/(4a^) < 1), there is no upper critical temperature, 
and the phase separation takes place even at infinitely high temperatures 
(Tc — >• oo). If the 'microscopic' end cap and junction energies, eg and 



themselves depend on temperature or density (e.g., in microemulsions WW )i 
the phase diagram may become more complicated. Note that the phase tran- 
sition described by Eq.(p7D shows reentrant behavior. At high temperatures, 
there is no phase separation, as usual. At very low temperatures the param- 
eter Kq = e~''^l'^ tends to zero and there is no separation either^ as shown in 
Fig. m, because the number of thermally generated junctions is too small to 
drive a phase separation. 

The phase transition discussed here is of purely entropic origin, which 
augments a qualitative claim by De Gennes that the cross-links are ther- 
modynamically equivalent to attractive interactions The origin of this 
phase transition lies in the fact that although the translational entropy of 
the chains is lower in the dense phase, the total entropy is still higher due to 
high entropy of the self-assembled junctions, which are abundant in the dense 
phase. It is important to emphasize that the presence of junctions induces 
an attraction between all monomers and not only between the junctions, as 
we shall show in Sec. ^. Also, the phase separation discussed here is not the 
same as the percolation (connectivity) transition, where an infinite cluster 
appears, although at both transitions a macroscopic connected network is 



formed. As we shall see in Sec. p.2| , the percolation transition, at which an 
infinite branched cluster appears in the system, is not a thermodynamic, but 
rather a structural transition. In some cases, the percolation transition can 
be masked by the end-junction phase separation. It is also interesting to note 
that for the values of and ej at which the phase separation transition is 



absent (/iq > Kq), there is also no maximum in as follows from Eq . (p6| 

3.2.1 Osmotic pressure and the phase coexistence line 

The osmotic pressure 11 can be calculated from the free energy / of Eq. 



Note that the coefficient before /ig-dependent (ends) term is positive due 
to the fact that the exponent in (j)^^"^ is less than one, while the sign of 
the junctions term is negative. This emphasizes that fact that, although 
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both ends and junctions terms are negative in the free energy, the junctions 
term is thermodynamically equivalent to an attraction, while the ends term 
produces an effective repulsion. The location of the phase coexistence line is 
determined as usual requiring that the chemical potentials A* = f| and the 
osmotic pressures 11 in both phases are equal: 



Numerical solution of these equations is shown in Fig. |l|. 

3.3 Four-fold and higher functionaUty junctions: gela- 
tion vs. vulcanization 

The analysis of the previous section can be generalized to systems with cross- 
links of arbitrary functionality z, by the addition to the Hamiltonian H3 of 
Eq.(|l^) a term that couples z spins. 



with Kz = J^/^e"*^^ (cf. Sec. ; being the energy of the z-fold junction. In 
this case the resulting Helmholtz free energy per unit volume then takes the 
form 



where ai and are two numerical, lattice dependent prefactors reflecting 
the local ends and junctions configurations, respectively. There are two 
interesting points to be noted. 

(z) There is no thermodynamic phase separation for z > 4, as can be 
readily seen from Eqs.(p7|,p9D , because the resulting free energy of Eq. ( pQ]) 
is always convex > 0). The attraction induced by the junctions weakens 
with increasing z. The topological percolation transition is, however, still 
present, as discussed in Sec. ^]2| 

{ii) When 2; = 4, the four- fold junctions contribution to the free energy is 
indistinguishable from the two-body attractions between the monomers (at 
least at the mean- field level), both contributing term ~ 0^ to the free energy 
for 0^1. The free energy is 



M0i) = /i(02) n(0i) = n(02) 




z times 



//T ^ (1 - 0)(ln(l - 0) - 1) - «i/io0'/' - «.e-^^0^/' 



(29) 



//T-(- 




1 1 



(30) 



e 
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where the effective junction energy e'^ has been redefined to absorb the 
internal degrees of freedom of a junction, expressed in the prefactor a of 
Eq.(^).The term ^e""^* corresponds to the x parameter of the Flory theory. 
At the point where €4 = 0, the excluded volume repulsion is balanced by the 
junction-induced attraction, a situation analogous to the 0-point of polymer 
solutions. The location of the compensation point is intuitively very clear: if 
€4 = 0, there is no cost for creating a junction and the chains can intersect 
each other at will (the fact that the junctions are four-fold is crucial!), thereby 
acting as completely Gaussian phantom chains. 

The three- and four-fold junctions correspond to two physically distinct 
situations. Three- fold junctions describe the process of gelation when the 
monomers simultaneously cross-link and polymerize in the reaction bath. 
The clusters formed this way exhibit mostly three- fold junctions, simply due 
to the fact that the collision of three monomers is more probable that the 
four- fold collision. In addition, in some systems (e.g. microemulsions [^ ) 
it can be shown that the energy of four-fold junctions is higher than that of 
three-fold junctions. The four-fold cross-links correspond to the process of 
vulcanization, when existing long polymer chains are cross-linked by addition 
of cross-link molecules, by irradiation or other means 0,^,|23|. Formation of 



three fold junctions is structurally inhibited by the nature of the cross-link 
which cross-links two preexisting chains. Contrary to gelation, these systems 
show only four-fold junctions that result from an intersecting pair of chains. 



4 Spatial variations and correlations. 

The advantage of the 'n = 0' model is that it can be used to predict the 
spatial correlations of the monomers, ends and junctions, in addition to the 
thermodynamic properties discussed above. These correlations are obtained 
from calculations of the second moment of the relevant probability distri- 
butions, such as the two-point correlations of the physical parameters, e.g., 
the density. The density-density correlation function can be measured in 
scattering experiments. Apart from being interesting per se, the correlations 
can also provide information about structural, non-thermodynamic transi- 
tions which one might expect in such structurally complex systems. Here, 
we generalize the mean field treatment of the previous section to include the 
possibility of spatial variations; this formalism allows us to calculate the cor- 
relations in the system. Each spin is approximated by its ensemble average. 
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Figure 3: Schematic illustration of the formation of a connected network 
form disconnected clusters upon increase in the monomer density 0. 

but a potential spatial variation in the value of the average spin (due for 
example, to a spatially varying external field) is allowed for (cf. Refs.@, 0). 
In this case one must allow for spatial variation of the constants, which are 
denoted now as Jj, hi, Ki each one labelled by its local spatial index. Ne- 
glecting the fluctuations of each spin about its local average value (S'j), as in 
Sec. H, we have 

S.S, = {S. - {S,)){S, - (S,)) - {S.){S,) + {S,)S, + {S.)S, ^ -{S.){S,) + {S,)S. + 
^ ■' 



where (Si) is a local spin average. Similarly we can approximate the triplet 
term as: 

S,S,Sk - -2{S,){S,){Sk) + {Sk){S,)S, + {S,){Sk)S, + {S.,){S,)Sk) 

Using the same procedure as in Sec. |^ and writing the average of the local 
spin as {Si) = Si, we find that the grand-canonical potential per unit volume. 
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uj{si} (cf. Appendix for details) is: 

1 1 1 ^""^ 1 

uj{si}V = 2 XI '^'■^^^^ + ^ XI -^^i^i^jSk - X ^^(^ + 2 '-^^ X "^^■^J + ^ X 3"^*"^ 

ij ijk i nrii j,k 



(31) 

where ^ is an unconstrained sum over all possible pairs of sites j, p, belong- 
ing to the same triple k. (cf. Appendix |8^ ). 
duj Xi 



dhi 1 + 



duj 



77 ^ ^ JiSi-Sj + —hiSi + — ^ ^ KiSi-SjSk = —SiXi (32) 



51n Ji 2 ^ 2 2 



From this we find: 



Xi 



Ji 



and Si = ^20.(1-0.) (33) 



(^) 



where Sj are sums over nearest neighbors of the site i. From these for- 
mulae, the Helmholtz free energy = c<jV{0j} + ^j0jln Jj can be cal- 
culated. 

For Kq = ho = (corresponding to infinitely long chains) the free energy 

is 



1/2 



(34) 

which we explicitly write here to emphasize the very non-trivial coupling 
between the different 0j's coming from the underlined term. For a spatially 
uniform system, 0^ = (p, and the F{(j)^} has only the single term (1 — 0) ln(l — 
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0) of Eq.(P^) because the terms proportional to vanish in the hmit of 
Kq = ho = (infinitely long chains) considered here for simplicity. 

The focus of the present discussion is different from the treatment of Ref. 



||23|| , which was restricted to a study of spin-spin correlations in the context 
of the percolation transition. In this section, we focus on the physically mea- 
surable density- density correlations and the thermodynamic properties. The 
location of the percolation line can also obtained within slightly modified for- 
malism, at least at the mean field level, allowing to relate the thermodynamic 
properties of the system to its structure. 

4.1 Density fluctuations 

Extending the mean-field theory to treat small, local density fluctuations, 
the free energy as a function of the local densities, F {(f)^} is systematically 
expanded to quadratic order in the density fluctuations related to the 
local average value of the density and the mean density by the relation: 
(j)- = + 5(j)i. 

F{0J = + J] In J, + J] S;,'6<f>M^ + 0{6<f>') (35) 



where Fmf(0) is just the mean field free energy of Sec. ^ Fmf(0) = ^/mf(0), 
S^j} = %^g^'^ is a matrix of the second derivatives of and we have 

used the fact that ^ = /i^ = In Jj. Expansion in Eq. (|35D is a completely 

general one, and the particular expressions for S^^} can be obtained from 
Eqs. (^,^). By the definition of the constant Jj, as an exponent of a 
local chemical potential, the density- density fluctuation function is given by 
{4'i4'k) ~ 0^ = mn J ainjfc which turns out to be the inverse of the matrix S~i}: 



" VdXnJk V 

where we have used the fact that ^^'^i^p^ = 0r A cumbersome but straight- 
forward calculation using the Eq. (^) (cf . Appendix |8]^ ) gives 
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^(5,,, + (1-20)^ M (36) 



]=nni 



93/^(2^ ./. (1-0) '2^" - <1 - 2*'^ 5: + - E ^.^) 

' j=nni j=nnni 

Writing the Fourier transform 



we get 

= J2S'\p}5<l>ip}5<l>i-p} and = 

p 

with S'(p) = ^ S'ifce*^*^'''"^*-'. From Eg . (|36D we find (for the simple cubic 

ik 

lattice) 

S-\p) = [A + B{cos{p,) + cos(p^) + cos(p,)) + (37) 



+C(cos(p^) cos{py) + cos(p2) cos(py) + cos(p^) cos{p^))] 



where 



- - '^^-)^(2^-^'^^)"^l 
C . ,2,-l)l^,2,)V^ 

and we have used the fact that the Fourier transform (FT) of the delta 
function is equal to unity and 
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where is the unit vector that points in the direction of the site j. On a 
simple cubic lattice, Pa is the component of the vector p in the direction a. 
As one can easily see from Eqs.(^, 0), the inverse of the structure factor at 
the zero wavevector 5 (0) = / S{r)dr is equal to the second derivative of the 
mean field free energy: 

^ ' dcp" 1 - ^ (20)3/V/2 (20)i/2g3/2 ^-^^^ 

This is an expression of the fluctuation-dissipation theorem in the grand- 
canonical ensemble. The variance of the mean number of monomers in the 
grand canonical ensemble is 

{N^) - {Nf = j dr{{(P{r) - 0(0))') = T\/0\t = ^(0) (39) 

where kt is the isothermal compressibility]^. Noting that 

g'FMF(0) _ ^ _ }_9li _ ^ _i 

902 ~ d(l)~ (f dv~ 03^^ 

(n is the osmotic pressure, w = 1/0 and yU = dFMF{(p)/d4> is the chemi- 
cal potential), Eq.(p8D becomes identical to Eq.(^ . Thus, the calculated 
density-density correlation function satisfies the expected thermodynamic 
sum rules. 

The structure factor S (p) = {S(j){p)5^{—p)) is experimentally measurable 
by scattering experiments (neutron, light. X-ray etc.) where the intensity of 
scattered radiation at the wavevector p is proportional to S {p) . If 5(0) is 
determined experimentally as a function of 0, the free energy -Fmf(0) can be 
found by the integration of Eq. (^^ . Obviously, Eq. ( |37D is not valid for the 
values of parameters at which the phase separation is observed, that is, in 

d'^ Fid)) 

the region where < 0. Rather, each macroscopic phase in a two-phase 

equilibrium has its own value of S{p), with the value of the monomer density 
appropriate for each phase. 

The expression (|37|) for S{p) is quite complicated, and its behavior is 
rather different for low and high densities. We next consider the nature 
of the predicted structure factor of Eq.(|37|) in the limits of high and low 
monomer density, 0. 
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4.1.1 Low density limit 



For small densities (j) the coefficients B and C are negative, and one can 
expand Eg . (1371) in a power series in the wavevector p. One finds that the 
coefficient of is positive so that this expansion is reasonable for wavevectors 
that are small compared with the inverse of the monomer size. 

S-\p) ~ 5-1(0) + a/ so that ^^^^ 



1 + aS(0)p2 

92Fmf(0),_i ^ (1-20) ,^ 3iro(«V3)15,^^,i/2 k 



5(0) = C ,^"^'"^ )-i and a = (1 + 7/^^ ^(20)1/^ _ 



20(1-0)' g3/2 4 ' (20g) 

Note that a is a non-negative, because ^^^^"^/a = ^) the half of the ratio of 
the number of ends to the number of junctions. In the present approximation 
of sparse junctions and ends,^ <^ 1 and a is non-negative. In the case of 

linear polymers without junctions, ^ is equal to the inverse average length 
of the chains, ^ = jf, which leads to the known correction to the scattering 
structure factor 0. 

In real space, the density-density correlation function is given by Fourier 
transform of the 5* (p). It follows from Eq.(^) that in real space ((50(f)50(O)) 
has the usual Ornstein-Zernicke form 



(50(r)50(O)) = J2 S{p)e''^"' ^ -e 



r 



with the correlation length ^ given by 



e = (41) 



The correlation length of density fiuctuations, ^, diverges at the spinodal line 
of the first-order junction-ends transition studied in the Sec.|3.2|, at which 



^ ^d^^^^^ ~ ^' ^^i^ expected for any first-order transition. Similarly, from 
Eq.(PD we see that the density fiuctuations at zero wavevector, S{p = 0), are 
also divergent at the spinodal, and, in particular, at the critical point. The 
divergence of the scattering intensity at the critical point should be observable 
in scattering experiments in the same way as the usual critical opalescence. 
We note in passing that in the absence of junctions, the structure factor, as 
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given by Eq. (HO) , reduces to the classical RPA result for polymer solutions 



It is interesting to note that the overall structure of a solution of self- 
assembled branched aggregates is not strictly self-similar, as we can see from 
Eq. ( ^0|) . There is at least one characteristic length in the problem, namely, 
the correlation length ^. At distances larger than ^ the correlations decay 
exponentially while for self-similar structures one expects algebraic decay of 
correlations. However, each branched cluster is expected to be self-similar. 



in analogy to percolation clusters |]2^. Note that both the correlation length 
^ and the scattering intensity at zero wavevector S (0) are both proportional 
to the second derivative of the free energy, This means that both 

of them are non-monotonic functions of the density, but have a maximum 
around the line 0g ~ 0^-, as follows from Eqs.(^ (pl), by taking the deriva- 
tives ||, The maxima of 5* (0) and of ^ as a function of are determined 
by the condition = ^^(^^"g^^ = 0, which is precisely the same con- 
dition which determines the location of the critical point (^~9^^l<^c ~ 
In particular, this leads to the conclusion that the maximum scattering in 
the region above the phase separation should be observed along a line which 
starts at the critical point and is determined by the condition ^ ^g^i^'^^ ~ 
( equivalent for the small densities to the condition 0g ~ 0^) , as shown in 
Fig. ^. This is in agreement with observations from scattering experiments 
from solutions of branched wormlike micelles [^(cf. also Discussion). 

It is instructive to compare ^ with another length scale that is present in 
the problem, namely, I, the mean distance between the defects (i.e. ends or 
junctions). As shown in Sec.^, 

7 20 

30, + 0, - 37^^0=^/2 + /.'o0^/2 ^ ^ 

and /(0) has a maximum around ~ h'^/ K'q. The behavior of / is quite differ- 
ent from that of ^, indicating that these two lengths are physically unrelated. 
The behavior of I is non-singular, while the correlation length ^ diverges at 
the spinodal line of the junction-ends transition. This emphasizes the fact 
that the junction-induced attraction has an effect on all the monomers in the 
system, as reflected in the behavior of ^; the effect of junction induced at- 
traction is not limited to the behavior of the junctions alone; this is reflected 
in the non-singular nature of /. The physical meaning of ^ is similar to the 
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'blob' size of the semidilute polymer solutions. Namely at distances smaller 
than ^, the behavior of an aggregate is that of a single self-avoiding branched 
polymer. 



4.1.2 High densities 



When (f) > 1/2, the coefficient (20 — 1) in Eg. (p^ becomes positive, which 
results in a negative coefficient in front of the term proportional to (in 
the expansion of the structure factor for small wavevectors) . Thus, the small 
wavevector expansion becomes meaningless, since it would give a negative 
structure factor for high enough values of p. Therefore, in this region one 
cannot use the expansion and one must retain the full form of the S{p). The 
exact location of the line where the behavior changes, = 1/2, is an artifact 
of the lattice model and should not be taken literally. Inspection of Eq.(|37D 
reveals that for small enough junctions densities, ( small Boltzmann factor 
Kq = e"^-'/^), for which 

{K'q = Koa/q'^^'^), the structure factor S (p) is a monotonically increasing 
function of the wavevector, p (for concreteness we consider the [1,1,1] direc- 
tion) for all > 1/2. 

This behavior of the structure factor at small junctions densities (weakly 
branched chains) agrees with the previous studies of the density fluctuations 
in concentrated polymer solutions and melts . Indeed, it is known that the 
structure factor of concentrated polymer solutions and melts behaves quali- 
tatively like that of simple liquids||3^. In particular, for small wavevectors p, 
the S (p) is an increasing function of p. This is related to the low compress- 
ibility of dense liquids. As discussed in the previous section, as the monomer 
density tends to one, the compressibility tends to zero. The magnitude 
of the density fluctuations, related to the compressibility by the fluctuation- 
dissipation theorem tends to zero as well 0. When the monomer density 
is high, but lower than one, the long wavelength (small p) density fluctua- 
tions are small, because large collective rearrangements of the monomers are 
structurally inhibited. However, short wavelength, local, rearrangements of 
the molecules are still possible. Combined together, these two effects produce 
a structure factor that increases with increasing wavevector, p. 
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In real systems, for the values of p larger than the inverse molecular size, 
one observes oscillations of the structure factor as a function of with the 
peaks corresponding to the first-shell, second-shell, etc., correlations of the 
positions of the adjacent monomers. Although the lattice model in the mean 
field approximation cannot be used to predict precisely the extremely low 
short scale features, such as the first shell peak observed in scattering exper- 
iments on polymer solutions, it adequately treats the qualitative behavior of 
the structure factor for wavevectors smaller than the inverse of the lattice 
constant. 

For higher values of Kq, K' > K*, (i.e. for higher junction densities) a 
peak emerges in S (p) { see Fig. ^) at a value of wavevector given by: 

Po = arccos [ - —J (43) 

indicating the presence of medium range correlations. This peak is quite 
different form the first-shell peak of simple fluids, which reflects excluded 
volume interactions on the molecular scale. On the contrary, the peak of 
Eq. (^^ reflects medium range structural correlations. By medium range we 
mean that the correlation range is larger than the lattice constant. The 
location of the peak depends on the values of the parameters, as discussed 
below, and it is not a lattice artifact. Because the peak appears only deep in 
the network phase, where the number of ends is negligible compared to the 
number of the junctions, one can neglect the /iQ-dependent term in Eq. (^7|) 
as a first approximation and write: 

indicating that the location of the peak moves towards lower values of p as 
K'q (20)^''^ increases. Fig. ^ shows the structure factor for two values of K'q-. 
below and above K*. Although Eq.(^4D was obtained for the [1,1,1] direction 
(Px = Py = Px), a similar result would be obtained for the isotropically 
averaged structure factor S (p) = J cPpS {p). The peak in 5* {p) should appear 
whenever there is a peak in any 'crystallographic' direction. In isotropic 
systems, the averaged structure factor S{p) is, of course, the experimentally 
measurable quantity, free of lattice and model artifacts. The predicted peak 
in the structure factor due to the presence of the junctions, is not a lattice 
artefact, because the location of this peak location is not related to the lattice 
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Figure 4: Structure factor as a function of wavevector p in the high density 
region, > i. Dashed hne shows the structure factor at low tempertaures, 
where the number of thermally-generated junctions is small. For higher tem- 
peratures, a peak developes in the structure factor, shown by the thick line. 

constant, but determined solely by the numbers of junctions present in the 
system. 

The precise size of the junction is arbitrary to a certain extent and can 
be adjusted to reflect the real physical features of any given system. We have 
defined a junction as a connection between three adjacent sites. However, 
one can also consider a larger junction, connecting next nearest neighbors, 
as in Fig. ^(b). Intuitively, one expects a more pronounced peak for larger 
junctions, because it amounts to an increase of the junction volume. Exact 
numerical calculations show that it is indeed the case. 

In the context of microemulsions, this peak can be related to the experi- 
mentally observed peak in the scattering from bicontinuous microemulsions 
(cf. Discussion). It is important to reahze that the absence of a peak 
does not, in principle, imply the absence of a network, (cf. Fig. ^). 
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5 Spin-spin correlations and topological struc- 
ture. 



So far we have been concerned with thermodynamic properties of a sys- 
tem and related phenomena such as equihbrium density fluctuations. They 



can be adequately described by a model presented in Sec. which maps 
the solution of equilibrium branched clusters onto Heisenberg model with 
anisotropic Si^iSj^iSk^i three-spin term. However, no structural information, 
e.g., concerning the formation of a continuous network, can be extracted 
from the model as it is formulated in Sec. |2.2|. Unfortunately, in order to be 



able to extract structural information about the network topology from the 
model, while retaining the exact correspondence between the spin model and 
the physical system of branched clusters, one has to resort to tensor order 
parameter, which results in a formally complicated theory |^3|]. Instead, we 
shall use a simpler, but less rigorous theory, applicable only in the mean field 
approximation, employed in this paper. 

Although the spin S has 'zero' components, one can still think about a 
'parallel' component Si pointing in the direction of the field h and a 'trans- 
verse' component S±. The rigorous mathematical procedure that justifies this 



is shown in Appendix. and amounts to calculating all the quantities of 
interest for a finite number of components n, and taking the limit n — 
subsequently. One can then devise a slightly more generalized version of the 
Hamiltonian |lO|, where the three-spin term takes the following form: 



ijk 



(45) 



where the sum is over distinct triples ijk and /5 is a numerical parameter. One 
can also include coupling between the 'transverse' components of the spin 
S, consistent with the symmetries of the system. This modification of the 
Hamiltonian has no effect on the thermodynamic properties and equilibrium 
fluctuations in the mean field approximation, discussed so far. Calculation 
of either the partition function or the two-point correlation function has 
contributions from the terms proportional to (3 that have at least one power 
of the average of the transverse component {S±). Since {S±) = 0, terms 
proportional to /3 do not contribute to thermodynamic quantities at the level 
of the mean field approximation. For (3 = the Hamiltonian (E^) reduces 
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to the previously used one. For /3 7^ the correspondence between the spin 
model and the equilibrium branched clusters is not exact. The clusters 
containing internal loops enter into expansion with the wrong weight due to 
summation over n — 1 'transverse' components. However, with the choice 
/3 = 3 all loopless tree-hke clusters are counted correctly (cf. Appendix 
p.4|) . Mean field approximation neglects the presence of intra-cluster loops, 
which formally correspond to correlations between different spins. This can 
be seen from the fact that the mean field free energy does not depend on 
the value of P in Eq. (^5|) . That is, it is independent of the precise way 
in which weights are assigned to the internal loops, which means that the 
intra cluster loops are neglected in the mean field approximation (cf. also 
Appendix |8.4| ). Thus, at the mean field level, the Eq.(^) with f3 = 3 can 



be used to calculate the transverse spin correlations which are related to the 
cluster size and percolation threshold. Therefore, we use the (3 = 3 model 
which neglects inter-cluster long-range correlations, responsible for the loops, 
consistently with the mean field approximation. 

Although the spin S itself has no physical meaning, it enters the calcu- 
lation of the various physical quantities, as we have seen in Sec. ^. Anal- 
ogously, the spin-spin correlation functions (Si^aSj^a), although unphysical 
themselves, can be used to derive physically relevant quantities such as the 
density-density correlation function. Another interesting, physical observ- 
able (e.g., by scattering experiments conducted in solutions of end-labelled 
chains) is the correlation between the end-points of the aggregates, which is 
related to the spin-spin correlation function by (cf. Appendix. p.4| ) 



{<Pe{r^)<PeirJ)) - <PI = = 5., + h[{S,S,) - {S,){S,)] (46) 

This can be understood geometrically by noting that the graphs which enter 
into (SiSj) are those with ends present both at the site i and site j.(cf. Fig. 
^ and Appendix |8.4| ) . 

Similarly, although S± itself does not have a direct physical interpretation, 
its correlator {Si^±Sk^±) does. If one considers the expansion of the correla- 
tion function {Si^±Sk,±) (cf. Appendix |8.4|) in powers of J, /i ( in the absence 
of junctions, i.e., K = 0) one can see that the only non-zero contributions 
arise from configurations that contain a single chain which starts at the site 
i and ends at the site k. Consequently, {Si^±Sk^±) measures the correlations 
between the ends of a single chain. As a matter of fact, even when K ^ 0, 
{Si^±Sk,±) is non-zero only when i and k are the ends belonging to the same 
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cluster (cf. Appendix and Fig. Thus, in general, the 'transverse' 
correlation function measures the correlations between chain ends belong- 
ing to the same aggregate. Similarly the 'longitudinal' correlation function 
{Si,iSk,i) measures the end-end correlations between any two ends, even if 
they belong to different chains or aggregates. The density-density correla- 
tions can be calculated from the 'longitudinal' correlator {Si^iSk^i) using Eq. 
(|3^ ) that relates (p^ to Si. This calculation yields result identical to those 
obtained in Sec. |873| . 

In this section, we present the results of a calculation of the spin-spin 
correlators (Si^aSj^a) using the same local mean field approximation which 
was used in Sec. ^ to find the density- density correlation function. Each 
spin is approximated by its local ensemble average: Si^a = (Si^a) local (with 
a = {1,_L}), which can, however, be different for spins at different sites 
labelled by i. Eventually, we will be interested in the deviations 6si of the 
local average spins from the spatially averaged value 5*. The full calculation 
is presented in Appendix. |8.4| . The spin-spin correlator is obtained from the 
partition function Z using the relations: 



d'^ ds ■ 

Using the results of Eq. (^) , a tedious but straightforward calculation yields 
the Fourier transform of the spin-spin correlator, C5 = ^^e*^'-'"'"^''-* (Ssi^aSsj^a) 

Aa 

fO. _ /AY) 

^ l-A'^{Jg,{p) + &Ksg2m ^ ^ 



where (cf. Appendix. ^]4| 



and where 



A« = (l-20)(l-0) 
A(^) = (l-d)) 



unit cell second cell 



gi{v) = E ^ 



unit cell 
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Figure 5: At low densities, < ^, the structure is a decreasing function 
of thewavevector. The scattering intensity at zero wavevector and the cor- 
relation length diverge at the phase stability coundary of the junctons-ends 
transition, shown by the thick closed line. In the region above the phase 
separation, the scattering intensity has a maximum along the dashed line, 
starting at the critical point. For high densities and low numbers of junc- 
tions, the structure factor is an increasing function of the wavevector, while 
for higher numbers of junctions present, a peak develops in the structure 
factor 

Consequently, the spin-spin correlation function in Fourier space is: 



In the following we show these results can be used to determine various 
physical characteristics of the system. 

5.1 Longitudinal correlations and the defect correla- 
tion function 




39 



where S (p) is the same as in Eq. (P7| ) for density- density correlations. Apart 
from the prefactor and the constant ( equal to 1) that ensures that the corre- 
lation of an end with itself is not counted, the dependence of (50e(p)^0e(~P)) 
on the wavevector p is identical to that of the full density-density correlation 



function of Sec. |8.3| . This is not accidental, but rather is a direct conse- 
quence of the fact that the density-density correlation function is related 
to spin-spin correlation function via Eq.(^). In particular, the correlation 
length at small densities ^ is the same as the correlation length of the total 
density fluctuations. The same is true of the junction-junction correlation 
function 

(W^Wj)) - ^3 = = '^^^■^"^ +6^^ ^dlC 

which can be calculated in a manner similar to that used to calculate (SiSj), 
and the subscript (K) indicates that the summation is over all possible cou- 
ples of the spins belonging to the original triple i, k, m. 

The fact that the correlation length of total density fluctuations equals 
that of the ends/junctions fluctuations emphasizes the point that the junc- 
tions induce an overall effective attraction between the monomers and not 
only between the defects (i.e., ends and junctions). In particular, the cor- 
relation length of junction-j unction correlations is the same as that of the 
monomer density-density correlations. The correlation length, related to the 
thermodynamics of the system, is unrelated to the mean distance I between 
the junctions which reflects the structural and topological properties. 



5.2 'Transverse' correlations: the structure of a single 
aggregate and percolation 



As have been mentioned, the 'transverse' correlation function measures the 
autocorrelations (denoted by the subscript auto) between the ends of the 
same cluster 

Because A-^ = (1 — 0) is always positive, we can safely expand the expression 
for the spin-spin correlator of Eq.(|47D up to second order in the wavevector p. 
Substituting the mean field values of J and K from Eq. (|2T|) , one gets 
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Figure 6: 'Longitudinal' spin-spin correlator measures the correlations be- 
tween any two chain ends. The 'transverse' correlator measures the correla- 
tion between the ends of the same chain. 



^ l-{l-(f))A^{Jgi{p) + 6Ksg2{p)) ho/ {2(f)y/^ - 3Koa{2(j)y/^ + bp^ 



where 6 is a non-singular function of the parameters. In other words 



with L = (20/((/), - 3(/)j))V2. This means that the correlations between the 
ends belonging to the same aggregate are negligible for r » L. For large 
separations r, the dominant contribution to (0g(r)0g(O)) comes from the 
largest cluster. In other words, the average size of the largest cluster is 

R^{^lt\"' (49) 
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When (f)j = the mean size of the largest diverges, and an infinite cluster 
is formed. This process is known as percolation and the emergence of an 
infinite network at (f)j = |0g agrees with other mean-field level studies of 
percolation p3[,[|]. Substituting for cj)^ and 0^ from Eq. ( p3D we obtain the 



following equation for the percolation line in the {4>,T) plane (cf. Fig. |l]): 

The length scale L = R has no thermodynamic meaning and purely refiects a 
topological (infinite network vs. disconnected clusters) property. At infinite 
temperature, the junction- induced interaction becomes irrelevant and 0p is 
the percolation threshold for the usual lattice percolation, studied in many 
classical works |2^. The value 0p = ^ (equal to 2{g-2) cubic lattice). 



obtained here is not inconsistent with the usual result 0^ = l/(g — 1) for the 
site percolation; the difference stems from the fact that in our model, there is 
percolation of two different objects: the monomers and the junctions, which 
do not occupy the sites of the same lattice; in our model, the monomers 
occupy the vertices of the lattice, while the junctions occupy the interstitial 
positions. 

Although the percolation transition, where an infinite cluster is formed 
is not related directly to the junctions-ends transition discussed in Sec. |3T^ , 
the junction induced interaction infiuences the shape of the percolation line 
at low temperatures, in agreement with earlier studies of the percolation in 



the interacting lattice gas [^. The exact value of a, and consequently of 0^, 
should not be taken too seriously, because the lattice approximation is not 
adequate for description of the system at the scale of an individual junction. 
Although, the calculation shows that the percolation line cuts the coexistence 
line of the phase transition on the left side of the critical point (cf. Fig. |I]), 
in real physical systems, the intersection point can be also to the right of 



the critical point |^|. Rigid systems can have effectively different a, as 
discussed in Sec.|^. Thus, depending on where the percolation line intersects 
the coexistence curve, in general, the coexisting phases discussed in Sec. 
^ can be either two network phases, two phases of disconnected, branched 
clusters, or a network coexisting with dilute chains. 

We now discuss the details of the cluster evolution as a function of 
monomer density (p. Recall that the mean distance between the defects (i.e. 
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junctions or ends) 

r 20 



+ 30, 



(50) 



is a non- monotonia function but has a maximum at 0g = 30,-; this is precisely 
the condition that determines the percolation threshold. Noting that Eq. (^9D 
imphes that the cluster size increases monotonically with increasing density 
0, we conclude that there are two competing processes in cluster formation: 
the lateral growth and the 'filling' in from inside. As the monomer density 
is increased, for small densities up to the percolation threshold, = 0p(T), 
the clusters become larger but more sparse. At 0p(T), an infinite cluster 
is formed and the process is reversed: clusters (including the infinite one) 
become denser, being predominantly filled in from inside. Although both 
below and above 0p(T) the lateral size of the clusters increases, the rate 
of lateral growth with increasing monomer densities lower than the rate of 
internal densification above (f)p{T). 



6 Rigid and semiflexible chains 

Solutions of rigid and semiflexible chains present two other broad classes of 
systems which may show self-assembly and formation of ends and junctions, 
in a manner similar to that of flexible chains as discussed up to now. As an 
example, one can think about solutions or gels of actin, a biopolymer forming 
the cytoskeleton and involved in cell locomotion. Although it is difficult to 
introduce bending rigidity into a lattice model such as ours, the qualitative 
effects can be easily understood, at least at the mean field level. For fiexible 
chains, as in our model, the lattice size is usually taken to be the persistence 
length (Kuhn segment for polymers), defining an effective 'monomer' size, 
while for semi-flexible chains, the lattice size would be the true size of molec- 
ular monomers, comprising the chains, the persistence length being much 
larger than the lattice size. On very general grounds, each defect (i.e. an end 
or a junction) contributes the energy —ksT to the free energy. Consequently, 
the defect-dependent part of the free energy in a system where non-conserved 
defects are present is 

F/T = -0defects = 
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in accord with Eq. P^ ) The difference between flexible chains and rigid 
ones lies in the dependence of 0g and 0^ on the total monomer density 0, 
which can be understood on the basis of the simple probabilistic argument 



outlined in Sec. |3.2| . Two ends can join to form a bond. In the case of 
flexible chains any two ends which are neighbors on the lattice can form a 
bond, irrespective of the orientation of the adjacent links (see Fig. 0). In 
the case of rigid rods, only rods that are co-linear, can coalesce. Therefore, 
taking into account that the coalescence of two ends lowers the energy by an 
amount 2ee, so that the Boltzmann factor for the probability of coalescence 
is proportional to e^^"", we find that in equilibrium 



c = q flexible chains 
c = 1 rigid rods 

The same argument applies for a formation of a junction from three ends 

a = g(g — 2)/3 flexible chains (cubic lattice) 
a = 1 rigid rods 

In general, for semiflexible chains 1 < c < g; a < g(g— 2)/3. In the continuum 
limit q ^ 4:71 for flexible chains. Therefore the most general form of the mean 
field free energy of Eq.(p^ can be written: 

F/T ~ (1 - 0) ln(l - 0) - (^)i/2g-e. _ ^rtf/2^-e, 

c c 

with c and a depending on rigidity of the chains and of the junctions. In 
general, the degree of the flexibility of the junctions is not necessarily equal 
to that of the chains themselves. The former may depend on the microscopic 
properties of the cross-linker molecules as in the case of gels or on the prop- 
erties of the surfactant molecules in the case of micelles and microemulsions, 
and other system-specific details. 

The effect of the rigidity on the correlations is more difficult to treat 
quantitatively in the present framework due to formal difficulty of introducing 
the bending rigidity into the model. In particular, the influence of branching 
on the nematic transition, present in the solutions of rigid rods remains to be 



investigated. An heuristic theory 41 1 might be more useful for this purpose. 
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7 Discussion 



We have presented a generic model of self-assembling chains which can branch 
and form networks with branching points (junctions) of arbitrary function- 
ality. The model rigorously maps the partition function of a solution of 
branched, self-assembling, mutually avoiding chains onto that of a Heisen- 
berg magnet in the mathematical limit of zero spin components. The model 
has been studied in the mean field approximation, which neglects the pres- 
ence of intracluster loops. It is found that despite the absence of any specific 
interaction between the chains, the presence of the junctions induces an 
effective attraction between the monomers, which in the case of three-fold 
junctions leads to a first order reentrant phase separation between a dilute 
phase consisting mainly of single chains, and a dense network. The reason for 
this entropic transition lies in the observation that although the translational 
entropy of the chains is lower in the dense phase, the total entropy is still 
higher due to the entropy of the large number of self-assembling junctions 
in the dense phase. The prerequisite for the phase separation is that the 
energy of a junction is lower than the energy a free end: ej < |ee, that is, the 
formation of junctions is energetically favorable (both ee, are considered 
positive and are measured relative to the energy of a monomer in the middle 
of a chain). The model then modified in order to study topological proper- 
ties of the system at the mean-field level. Independent of the junctions-ends 
transition, we predict the percolation (connectivity) transition at which an 
infinite network is formed. The percolation transition partially overlaps with 
the junctions-ends transition (cf. Fig. |1|). This transition takes place at the 
point where the density of the free ends is three times the junctions density, 
0g = 3(f) j, which at the same time is the point where the mean distance be- 
tween the defects (i.e., junctions or ends), T attains its maximum. This result 
agrees with other mean field level studies of percolation of tree-like clusters 
This result means that up to the percolation transition the branched 
clusters grow laterally, predominantly by addition of the new branches at 
the periphery of a cluster, while at the percolation transition, an infinite net- 
work is formed and the process is reversed: the clusters grow mainly due to 
the filling of the internal mesh. The percolation transition is a continuous, 
non thermodynamic transition that describes a change in the topology of the 
system but not a thermodynamic phase transition. Our treatment which 
predicts both the thermodynamic phase equilibria as well as the spatial cor- 
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relations in the system allows us to treat both the phase separation and the 
percolation threshold within the same framework, at least at the mean field 
level. 

The predicted density-density correlation function has a usual Ornstein- 
Zernicke form at low monomer densities. The correlation length of the density 
fiuctuations ^ is divergent at the spinodal line of the junctions-ends transi- 
tion and is a non-monotonic function of the monomer density. C, attains its 
maximum around the point where the number of ends equals the number of 
junctions, 0g ~ (pj. It is important to emphasize that it is ^ which is mea- 
sured in actual scattering experiments, and not other lengths of the problem, 
such as the mean cluster radius R or the mean distance between the defects 
/. The zero-wavevector scattering intensity, S (0) also has a maximum as a 
function of (f) along the same line, starting at the critical point. This effect 
is a direct consequence of the first-order phase separation between low and 
high density phases (cf. Fig. |^) 

We also predict the emergence of the medium range correlations at high 
monomer and junction densities, as refiected in the predicted peak in the 
structure factor, signifying the onset of structural correlations in the system 
(cf. Fig. D. 

The theory presented here has many physical realizations. We now discuss 
the implications of the generic results for several specific cases. 



7.1 Physical gels. 

Physical gels include a broad class of systems, consisting of long polymer 
chains that are reversibly cross-linked. The cross-links can break and reform 
under the infiuence of thermal fiuctuations. Examples of these materials are 
numerous, many of them of practical interest p6[. The chains that comprise 



the basis of the gel, can either be self-assembling or chemically bonded ||22 
With increasing monomer density, the gelation threshold is reached where 
a connected network (gel) is formed. Despite intensive research in the past 
two decades, several features of the physical gelation process are still under 
debate. One question concerns whether the gelation transition (i.e., the 
transition from disconnected, branched clusters to a macroscopic, connected 
network) is a structural transition or a thermodynamic one. Another question 
of interest is: if the gel transition is indeed a thermodynamic one, what is the 
order of the transition ^ . Our results show that actually there are two 
independent transitions. The first one is the classical percolation transition. 
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located at the point where the density of the cross-hnks is equal to one third 
of the density of the free ends in agreement with classical result due to Flory 
0]. An infinite connected cluster is formed at the percolation threshold. 
This transition is structural and is not connected to any thermodynamic 
singularity. 

For three-fold cross-links another, a thermodynamic transition is also 
present at small densities, and partially overlaps the percolation line (see Fig. 
IH). This transition is the first order phase transition due to inter- monomer 
attraction induced by the presence of the junctions and not due to the aggre- 
gation of stickers (cross-link molecules). Our model shows that the phase sep- 
aration can occur under good solvent conditions, and not only in B-solvents. 
An advantage of the present model is that it starts from a rigorous descrip- 
tion of a solution of self-avoiding branched chains and does not involve any ad 
hoc hypotheses about interactions between the monomers, cross-linkers and 
solvent molecules. An interesting problem is the difference between gelation, 
when molecular segments are allowed to polymerize and cross-link simulta- 
neously in the reaction bath, and vulcanization, when preexisting polymer 
chains are cross-linked by irradiation, addition of linker molecules, or other 
means. In addition to the fact that one expects large spatial inhomogeneities 



in the case of vulcanization we argue that an additional major differ- 
ence is the functionality of the cross-links. In the case of gelation, most 
of the junctions will be three-fold, simply due to the fact that the collision 
of three segments is more probable than a collision of four. In the case of 
vulcanization, on the contrary, the junctions will be four-fold, being formed 
by cross-linking two pre-existing chains. So, in certain cases gelation will be 
observed experimentally as a first-order thermodynamic transition, while vul- 
canization is always continuous non-thermodynamic transition. As a matter 
of fact, we have shown that the presence of four-fold junctions is thermo- 
dynamically equivalent to decreasing the quality of the solvent, in accord 
with the results obtained by other theoretical methods and Monte Carlo 



simulations [26 



The long-standing question about the presence of closed, intra-cluster 
loops in the pre-gel clusters and in the gel phase, becomes irrelevant in the 
present formulation. In our model, nothing prevents the intra-cluster loops 
from forming and their number is determined by the distributions of ends and 
junctions. However, we note that the question of the influence of the loops on 
the thermodynamic and structural properties cannot be addressed properly 
in the mean field approximation, used in the paper, because it effectively 
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disregards the loops. 

Experimentally, the gelation transition is commonly observed by mea- 
suring a viscosity increase due to the formation of the macroscopic net- 
work. Therefore, the nature of the experimentally observed gelation tran- 
sition will depend on the path in the phase space along which the transi- 
tion is approached. For example, decreasing the temperature at constant 
monomer density along the arrow b of Fig.|I| will be produce a continuous, 
non-thermodynamic gelation transition at the point where line b intersects 
the percolation line. On the contrary, a temperature decrease along the ar- 
row a will produce a first order phase separation between a phase of dilute 
chains and a connected network (gel). 



7.2 Microemulsions 

Microemulsions consist of domains of two immiscible fluids (typically, water 
and oil) stabilized by a surfactant 0,0. One theoretical approach to their 
understanding is to map the oil-in-water dispersion onto a lattice model in 
the following manner. When the spontaneous curvature is small, the 
relevant length scale is the persistence length of the surfactant film; this 
length is determined by the bending rigidity and the thermal fluctuations 



4^ . For systems where the spontaneous curvature dominates (this occurs 
when the spontaneous curvature is larger than the inverse of the persistence 
length), for a certain range of spontaneous curvatures, one finds that the oil 
or water domains are either spheres or elongated tubes of radius r ~ Cq ^. For 
oil internal systems, the radius of the oil domains can also be kept constant 
by maintaining fixed ratio of the surfactant density to that of oil (in 
a water internal system it is given by the surfactant to water ratio). The 
application of the general results of this paper to these cases requires that 
one understands the mapping between the microemulsion and our lattice 
model. Taking r as the lattice constant, and considering an oil internal 
system, the microemulsion can be mapped to a solution of self-assembling, 
branched chains (oil domains-^ 'monomers'; oil+surfactant fraction— (f)) as 
described by our model. Obviously, the analogy breaks down at oil volume 
fractions close to one, because our model does not allow the merging of a two 
elongated thin adjacent domains into a wider one. 

Using this analogy, the evolution of spontaneous curvature dominated 
microemulsions can be qualitatively described in the context of our model and 
the results presented above, as follows. At low oil fraction, the microemulsion 



48 



consists of spherical droplets of radius r ~ Cq ^. At a certain oil fraction a 
'polymerization transition' [^^] takes place and cylinders of radius r are 
formed which eventually branch. With increasing oil fraction at a certain 
oil fraction (j)p(T) ~const xe^^''~^^''^'^ , a percolating oil domain is formed, 
indicating the onset of the bicontinuous microemulsion (cf. Sec. [4.1.2|) . At 
low temperatures, this transition is masked by a first order, phase separation 
transition due to junction- induced attraction between the oil domains (Sec. 
p.2|) . This first-order transition is present only if the branching points are 
three-fold, which seems to be relevant to microemulsions. In addition to the 
fact that three-fold branches are more probable statistically, the energy of 
the three-fold junctions is lower than that of four- and higher fold junctions, 
as can be shown from calculations of the bending energy of the surfactant 
bilayers in the junction region [Q. 

The structure factor (for scattering from the bulk water or oil domains - 
i.e., bulk contrast experiments) has a simple Ornstein-Zernicke form in the 
droplet microemulsion region (which includes also elongated but disconnected 
cylindrical droplets) with a correlation length that diverges at the spinodal 
line of the first-order transition. In the bicontinuous region, the structure fac- 
tor is an increasing function of the wavevector, as expected in any high den- 



sity system (cf. Sec. [4. 1.2]) . With increase in the number of junction, either 
due to increase in the monomer density, or the temperature concentration 
a peak appears in the structure factor, indicating the emergence of medium 
range correlations between oil domains (cf. Fig. |^). The peak is predicted to 
appear whenever the condition (p > ^e"^'^^/"^ is satisfied. The peaked struc- 
ture factor is a distinctive characteristic of bicontinuous microemulsions, as 
borne out by many experimental studies It is important to empha- 

size, that our model predicts that the peak occurs only when (p > max[|, 0^] 
which means that, in principle there may be network microemulsions that do 
not show a peak in the structure factor ( the region between the percolation 
line and vertical = 1/2 line Fig. ^). The fact that the peak appears only 
for > 1/2 should not be taken too seriously. The exact value = | is most 
probably an artifact of a lattice construction, refiecting the monomer-hole 
symmetry. More important is the qualitative prediction that the peak in the 
structure factor is a consequence of the correlation between the oil domains, 
induced by the presence of the junctions. 
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7.3 Wormlike micelles 



Surfactant molecules in aqueous solutions can form long, cylindrical micelles 
in certain regions of the phase diagram. In the cylindrical phase, each micelle 
consists of a large numbers of surfactant molecules. The micelles are polydis- 
perse in size, the equilibrium distribution being determined by the interplay 
between the entropy (which favors small micelles) and the 'cap' energy, which 
is determined by the geometry of packing of surfactant molecules at the ends 
of a micelle. Due to the bending rigidity of the surfactant layer, the mi- 
celles are stiff up scales of order of the persistence length Ip |^ . On the other 
hand, for length scales much larger than /p, the micelles can be considered 
as flexible chains. Therefore, for long enough micelles (/io/a/20 -C 1), the 
lattice constant of the equivalent Heisenberg model can be taken to be equal 
to the persistence length Ip. In the opposite case of short chains, on can take 
a surfactant molecule as a 'monomer' size, with the reservations described in 
Sec. |. 

It has been suggested on the basis of rheological, conductivity and dielec- 
tric polarizability experiments, that the micelles can branch, and at a certain 
point the system transforms to a connected micelle network. The energy of 
the junctions and the ends in this system can be varied by changing the 
salt concentration. Theoretically, it was proposed that the ' n = 0' model 
might describe the formation of a network of wormlike micelles [pi] , ^ . Our 
results substantiate these suggestions and indeed predict the formation of a 
connected micellar network, which can partially overlap with the first order, 
phase separation phase transition. 

Results of neutron scattering performed on wormlike micelles solutions 
show that the correlation length, ^, of the density fluctuations is a non- 
monotonic function of the density in the proposed branched state I^S], which 
might be explained by the results of Sec. [4.1.1| . These experimental results 
were interpreted in Ref. in terms of the mean micelle size (analogous to / 
in our formulation). However, it is important to realize that, although both ^ 
and I exhibit a maximum as a function of the density, it is ^ which is actually 
measured in scattering experiments, while I conveys purely geometrical in- 
formation, not directly measured in scattering experiments. The correlation 
length ^ has a maximum around (f)^ = (pj which is a direct consequence of the 
junctions-ends transition, discussed in Sec. It should be noted that the 
maximum of ^ is not directly related to the network formation, which takes 
place at 0g = 30 -, as well as I. 
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7.4 Dipolar and magnetic fluids and colloids 



Dipolar fluids or colloids consist of dipolar particles, bearing electric or mag- 
netic dipole. For, example ferrofiuids. consist of magnetic (metallic) particles 
immersed in an inert fluid. The molecules of real liquids, such as acetone 
(CsHqO) or dimethylsulfoxide {C2H3SO) also bear electronic dipolar mo- 
ments. The interaction energy of two dipoles has a minimum when the 
dipole moments are coUinear |TI|. This means that dipolar particles have 



a tendency to aggregate into chains |46|, so that all the dipole moments of 



the particles in the chain are more or less colinear. However, it is important 



to realize that these chains can also branch The energy of the parti- 

cles at the branching point is obviously higher than in the middle of a chain 
and can be calculated form electrostatics (or magnetostatics in the case of 



magnetic dipoles) |^^. Disregarding the long-range part of the dipole-dipole 



interactions, the system of both finite ended and branched dipolar chains can 
be mapped to the model presented in this paper and our predictions for the 
phase separation, scattering and percolation would apply. This mapping is 
justified in dilute systems such as ferrofiuids. or magnetic colloids. In dense 
molecular liquids, it probably breaks down. In order to understand the prop- 
erties of molecular dipolar liquids, one may have to take into account specific 
interactions between the molecules, such as hydrogen bonding. 

7.5 Actin networks 

Actin is a self-assembling biopolymer, which supports the membranes of bi- 
ological cells. Sol-gel transition of actin solutions and networks, induced by 
various proteins, (e.g., myosin) plays an important role in cell movement 
and locomotion. The properties of actin networks in vivo have been stud- 
ied extensively by biochemical and molecular biology methods 0. Recently, 
physical properties of actin solutions and gels have been studied in vitro, 
in rheological and scattering experiments The experimentally obtained 
phase diagram is strikingly similar to Fig. |I|, where at the percolation line 
there is a transition from and entangled network to a connected 'microgel' 
state. However, it was observed that at low temperatures the dense phase 
at the coexisting curve consists of 'bundles', where several parallel actin fil- 



aments are tightly bound together by the cross-links [^ . 

Actin filaments are rigid, with the persistence length as high as ten mi- 
crons. However, the general results of this paper should apply, with the 
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reservations made in Sec. ^, which accounts for the similarity between the 
predicted phase diagram and the experimentally found one. The formation 
of the bundles is probably due to nematic interaction between rigid filaments 
and cannot be captured in the present model. However, one can estimate 
the condition for the bundles formation using a crude qualitative argument. 
For this purpose, the network phase can be viewed as a collection of rigid 
filaments of average length e^^^'^ (j)"^^"^ (cf. Sec. p.l.l]) . The orientational 
order appears in the solutions of rigid rods when > const/L, which means 
that for > e~^'^j/^ one might expect appearance of the bundles. 
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8 Appendices 



8.1 AppendixA: ' n = 0' model 

The formal aspects of the 'n = 0' model have been studied extensively 



r2||l9|. For the sake of completeness we review the derivation here. Let us 
calculate the generating function Tr e'^* for an n-component spin, s, on a 
d-dimensional lattice, where the averaging is over solid angle of the vector s, 
subject to condition J2a = ^- 

The cumulants of s*can be obtained by differentiating the generating function 
with respect to k. Now, supposing that k = {k, 0, 0) 
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The integrals in step two of the above equations, although formally divergent, 
should be treated as generalized functions. The resulting integrals are known 
Putting V = 1 — n/2, we find 



9n{k) 



n I -1 n 

7^2^-^77,2 



(fcv/^/2)t-i 



where Hi,^\x), h'P (x) are Hankel functions of the order p, and Ip{x) is a 
modified Bessel function of the order p. Now, up to second order 

Keeping in mind that in our case p = | — 1, it follows 



In particular. 
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(51) 



Due to the fact that the generating function has only quadratic terms in k, 
cumulants higher than second order vanish. This is because the higher order 
cumulants are related to higher order derivatives of the generating function 
with respect to k and these all vanish, as seen from Eq. (|5TD 



8.2 Appendix B: monomers, bonds and junctions on a 
lattice 

8.2.1 Relation between the densities 

The chains consist of monomers, occupying the vertices of the lattices, con- 
nected by bonds. Each junction is not occupied by a monomer but is formed 
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by three unphysical 'bonds', connecting three adjacent monomers as illus- 
trated in Fig. ^ The bonds forming the junction are not counted by the 
model. The monomers comprising a chain fall into two categories: internal 
ones, attached to two bonds, and the monomers at the ends of the clus- 
ter, attached only to one bond. The total number of bonds (including the 
unphysical ones) is therefore 

^^bonds = ^(2iV. + iVe + 3iV,) 

where Ni is the number of internal monomers, Nj is the number of junctions, 
and Ne is the number of ends; the prefactor ^ accounts for the double counting 
of each bond. The total number of monomers is 

N = Ni + Ne 

and the number of physical bonds, connecting a pair of monomers, is the total 
number of bonds less those unphysical bonds which are involved in junctions 

^bonds = ^bonds - ^^j 

Combining these equations together, one gets 

N = TVbonds + + ^N, 

8.2.2 Value of a 

Here we present the considerations which lead to particular values of the con- 
stants in mean-field equations of Sec. |^. In our model we consider each site 
of a lattice occupied either by a monomer or solvent molecule. Chain bonds 
are defined as the bond between two occupied neighbor sites. A junction can 
be formed when three ends are neighbors on a lattice as shown in Fig. 0. 
The sum in the three-spin term Yliijk-^^i^j^k in Eq- (|T8|), accounting for 
the junctions between the chains, is constrained over all distinct triplets of 
spins z, j, k. Tthe three-spin term can be rewritten as | KSi Yll^jk^ SjSk 
where Yli^^ unconstrained summation over all cites and Y2j k 

constrained sum over all possible pairs of sites j, k which together with i 
form a junction i, j, k. In the mean field approximation this term is fuirther 



54 



Figure 7: Two possible configurations of the triplet i,j,k: a) the j and k 
sites are the nearest neighbors of the site i; b) k is the nearest neghbor of 
the site i while j is the next nearest neighbor of i. 

transformed into ^'^i'iKSiY^^^\Sj){Sk) ■ On a hypercubic lattice, there 
are two possible configurations of the triplet i,j,k. The first one, is shown 
in Fig. 0(a): both j and k are the nearest neigbors of the site i. The second 
is realized when either j or /c is a next nearest neighbor of i, shown in Fig. 
0(b). A simple calculation shows that for a hypercubic lattice of a coordi- 
nation number q, the number of such pairs of the first kind is is q{q — 2) 
and the number of the pairs fo the second kind is \q{q — 2). Therefore, 
a = |(g(g — 2) + |g(g — 2)) = ^q{q — 2). For a hexagonal close-packed struc- 
ture or for the f.c.c. lattice, a similar calculation gives a = q (|l(a)). However, 
the exact numerical value of a is an artefact of the lattice construction and 
should be thought of as a phenomenological parameter refiecting the micro- 
scopic features of a junction, along with the junction energy ej. A particular 
value of a affects only the numerical values of the critical temperature and 
density but not the overall qualitative behavior. It can be adjusted to refiect 
physically allowable junction configurations. For example, three neigboring 
chains, coming perpendicular to the plane of a triplet i,j, k and connected 
by it, can hardly be interpreted as a junction in most of the experimental 
systems. For convenience, we have used a = |g(g — 2) in calculating the 
phase diagram and the structure factor. 



55 



Figure 8: Two possible alternative definitions of a junction, a) shows a 
junction on a hexagonal lattice, a = q/2; h) shows a larger junction, coupling 
the sites separated by two lattice constants. 



8.3 Appendix C: density-density correlations 

In this Appendix we show the derivation of the result (^) of Sec. ^ for the 
density-density correlation function. From Eq.(|33D 
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Now, 

Edsj 



J 



1 



Together with Eq.(^) this leads to Eq. (|36|) of Sec.^ 



8.4 Appendix D: spin-spin correlations 

This Appendix shows the detailed calculation of the spin-spin correlation 
function of Sec. |^. In principle, the scalar products of any two spins in Eq. 
( [T0| ) contain both longitudinal component. Si, parallel to the direction of the 
field /i, and a 'transverse' one, S±. For the two spin term it gives 

JSi ■ Sj = J (^Si^iSij + Si^± ■ Sj^±^ (53) 

and the three-spin term is chosen to be 

ijk ^ 

(54) 



The average value of the 'transverse' component, {S±), is equal to zero, be- 
cause the average magnetization is aligned along the field h. Consequently, 
it has no consequence for the thermodynamic properties, and was neglected 
in the calculation of the free energy in the Sec. ^. However, its fluctuations 
around its average of zero, are not zero, and convey physically important 
information. Neither Si nor S± has any direct physical meaning. However, 
physically relevant information can be extracted from the knowledge of their 
correlation functions. 
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Consider the correlation function (S'l.jS'ij) between the 'longitudinal' 
components of two spins located at sites i and j. As follows from the dis- 



cussion of the Sees. |^, p.2| the only terms entering its expansion Tr Si^iSij 
Q-H{Sk}/T g^j^g ^Yie same as of Eq. (|13]) with Si^iSij replacing the /i-terms 
hSi and hSj, This means that the correlator (Si^iSij) counts all the config- 
urations where there are chain ends at the point i and at the point j. In 
other words, {Si^iSij) measures the correlations between any two ends. 
Recalling that the ends density is 0e(rj) = hiSi^i, and subtracting the 5ij 
contribution corresponding to the correlation of an end with itself, we arrive 
at Eqs. (^61) (^8]), which can also be obtained directly from the relation 



(50, (r,)50, (r,)) 



dhj 



To calculate the 'transverse' correlation function, {Si^±Sj^±), where now 
S± signifies a component of S in one of the {n — 1) 'transverse' directions, 
one must focus on the second term of Eq.(p3D . Following the arguments of 



Sec. 2^, the only terms which enter into expansion of {Si^±Sj^±) in the limit 
n — >■ are of the form 

^i,±[JSi^±Sm,±JSk',±Sk,± [^KSm,±Sn,lSk,±]JSn,lSn',l--- (55) 

. ' 

junction 

...KSm',lSp^lSq^l JSj'^±Sj^±hSq^l]Sj^± 

that is, the clusters, consisting of the 'backbone' of the S±^s with branched 
'sidechains' of the S'l's, closed by Si^± and Sj^±. Inspection of possible terms 
of the form (^) shows Eq.(|53Dthat the sites i and j must belong to the same 
chain, because the Eqs. (|53| , |5^ couple either two transverse components or a 
longitudinal one with two transverse. In other words, the 'transverse' corre- 
lator measures the correlations between any two ends belonging to a same 
chain. It is seen from equation ( pS] ) that the junctions enter the expansion 
of the transverse correlator with the weight ^K. Since one wants to measure 
the correlations between the ends of a cluster belonging to the expansion of 
the original Hamiltonian ([T0|) , where all the junctions enter with the same 
weight K, the value (3 = 3 must be chosen. That is, the relevant part of the 
three-spin term in the Hamiltonian (p!OD has the form 

ijk 
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It is important to emphasize that the model given by Eq. (^) is suitable 
only for loopless clusters. Clusters containing loops enter the expansion with 
wrong weights, due to summation over n — 1 'transverse' components. How- 
ever, the clusters with loops are neglected in the mean field approximation 
and one can use the model given by Eq.(|5^). 



Performing the mean field approximation on these terms, as outlined in 



Sec. gives for the local spin averages : 

dlnZ xi 

Si.l = 



<9/i,,i l + ^[xl + xl 
dlnZ x± 
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where 

xi = hi^i + ^Jsj^i + 3{^^Ksj^iSk,i) (56) 

1 ^""^ 



3 

j,k 



with J2j is the sum over nearest neighbors of the site i and J2fk^ the sum 
over the pairs of sites j, k which together with i make the junction triplet 
i,j,k as in Sec. ^. The coupling between the parallel and transverse com- 



ponents of s in the three-spin term in Eq. (|56D has no consequence for the 
thermodynamic properties because the average value of 'transverse' compo- 
nent is zero: {S±) = s_l = 0. The terms containing products Si^±Sj^± have 
been omitted for the same reason. However, this coupling is important for 
correlations, as we shall see. We thus define the spin-spin correlation function 
Cl^\ whose 'longitudinal' component is given by: 



1 ^""^ 



,.2 



. ^ . d , X\ . . 1 X-^ 

^'^^^ = d^S + \[xl + xi] ^'^-=° = TTM ~ (1 + \xif 
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For transverse component, (bearing in mind that {S±) = s± = 0) we obtain: 



jGnrii 



Note that the couphng between s_l and Si in Eq.(^8|) does contribute to the 
transverse spin-spin correlation function. Substituting from Eq. (^) one finds 

= (i-20)(l-0) 
A(^) = (1-0) 

Recalhng that J^IH] = If Enn, + Ennn,' after Fourier transform C'l^'^ = 
^__gip(r- -r,)(-.W (a = {1, X}) the equations §^ and (H) yield the Eq. 



of Sec. It is important to emphasize that the model given by Eq. 
is suitable only for loopless clusters. Clusters containing loops enter the 
expansion with wrong weights, due to summation over n — 1 'transverse' 
components. In principle, one can include higher order terms, coupling dif- 
ferent 'transverse' components. However, as seen from the calculation above, 
in the mean field approximation, these terms do not contribute either to the 
free energy or to the spin-spin correlation function. That is, the mean field 
results are independent of the precise way in which the weights are assigned 
to clusters with loops. It follows that the mean field approximation disre- 
gards the loops, which contain information about long range intra-cluster 
correlations. Only the tree-like clusters give contribution to the mean field 
approximation. It is also possible to show by direct calculation from the 
rigorous Lubensky-Isaacson model that in the mean field approximation the 
number of loops vanishes. 



8.5 Appendix E: Exact model for percolation 

In order extract the information about the topological structure of the sys- 
tem, while retaining the exact correspondence between the expansion terms 
of the spin model, and the equilibrium branched clusters, one has to resort to 
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tensor order parameter or, equivalently, to use several coupled n = fields. 
The presentation in this section follows closely the derivation of Ref . with 



the major distinction that the junctions are not treated as point-like objects, 
which allows one to use a conventional normalization of the n = spins S. 
The Hamiltonian to be used is 



= ^ E E + E E ^"^T + E E ^M^".^"fe (59) 

i a=l (jj) a=l {i,j,k) ce=l 

where S'^'s are m n-component vectors with the normalization ^q,(5'")^ = 
mn where n —>■ 0; a goes form 1 to m. In this form, the Hamiltonian consists 



of m replicas of the Hamiltonian (|T3), and there is one to one correspondence 
between the terms in the expansion of Tr e~^^, and the equilibrium branched 
clusters. Due to presence of the additional parameter, m, each cluster in the 
expansion acquires an additional weight m. The partition function then reads 

Zt= J2 J^'K^' h^^m^^Af{Nj , N,, Nh, N,) 

Nj,N^,Nt,Np 

where Nj, N^, Nh, Nc are the number of junctions, ends, bonds, and clusters, 
respectively. The parameter m can be analytically continued to non-integer 
values and has the meaning of the fugacity conjugate to the total number 
of clusters in the system. For self-assembled branched clusters, the relevant 



value is m — i> 1, because the number of clusters is not constrained . 

The Hamiltonian (1) can be rewritten, making the rotation in the Tri- 
dimensional space of a-components, with the use of a set of m-dimensional 
orthogonal vectors vi; / = 0, m — 1, with Vo = ^(1, 1, 1) and Vi ■ Vu = 



Sw p3[. The new variables are then defined as: 

a * a 

= ^S^X for /^O 



Note that the normalization remains the same: Ylb=o^ (5*')^ = mn. The 
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Hamiltonian expressed in the new variables then reads 



m—1 m—1 
i 1=0 1=0 

m—1 



m , 



S^S^S'^ + I E + Perm(^, J, k)\+ E ^ii'i"SlS';si' 

1=0 / i^O,«VO,«'VO 



where aun" are numerical coefficients, depending on the value of m. It can 
be shown then that the correlation function {|S'°S'^) measures the correlation 
between any two endpoints and the (SlSj) measures the correlations between 
two endpoints belonging to the same cluster |]23[ . 

In the mean field approximation, the products of S'' with / 7^ do not 
contribute either to the free energy or to the two point correlation function, 
because only (S^) 7^ (they do, however, contribute to higher order correla- 
tion functions), and can be neglected. In particular, the fact that the mean 
field two-point correlation function is independent of the coefficients awi" 
(i.e. , of the precise way in which weights are assigned to loops), means that 
the mean field approximation takes into account only the loopless graphs 
contributing to the two point correlation function. With this in mind, the 
last term in Eg. (|60D vanishes, and in the limit m — 1 the Hamiltonian (2) 
becomes formally identical to the Hamiltonian (^5]) of Sec. ^ with S*' (/ 7^ 0) 
taking place of a 'transverse' component S^, 7 7^ 1- 
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